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Overview 

Precision measurements of electroweak quantities 
are carried out to test the Standard Model (SM). In 
particular, measurements of the top quark mass, rritop, 
when combined with precise measurements of the W 
mass, Mw^ and the weak mixing angle, sin^ 9w , make 
it possible to derive indirect constraints on the Higgs 
boson mass, Af/f , via top quark and Higgs boson elec- 
troweak radiative corrections to M^- Comparison of 
these constraints on Mh with the mass obtained from 
direct observation of the Higgs boson in future collider 
experiments will be an important test of the SM. 

In this report, the prospects for measuring the W 
parameters (mass and width) and the weak mixing 
angle in Run II are discussed, and a program for ex- 
tracting the probability distribution function of Mh is 
described. This is done in the form of three largely 
separate contributions. 

The first contribution describes in detail the strate- 
gies of measuring Mw and the W width, at 
hadron colliders, and discusses the statistical, theo- 
retical and detector specific uncertainties expected in 
Run II. The understanding of electroweak radiative 
corrections is crucial for precision W mass measure- 



ments. Recently, improved calculations of the elec- 
troweak radiative corrections to W and Z boson pro- 
duction in hadronic collisions became available. These 
calculations are summarized and preliminary results 
from converting the theoretical weighted Monte Carlo 
program into an event generator are described. The 
traditional method of extracting Mw from the line- 
shape of the transverse mass distribution has been the 
optimal technique for the extraction of Mw in the 
low luminosity environment of Run I. Other techniques 
may cancel some of the systematic and statistical un- 
certainties resulting in more precise measurements for 
the high luminosities expected in Run II. Measuring 
the W mass from fits of the transverse momentum dis- 
tributions of the W decay products and the ratio of the 
transverse masses of the W and Z bosons are investi- 
gated in some detail. Finally, the precision expected 
for the W mass in Run II is compared with that from 
current LEP II data, and the accuracy one might hope 
to achieve at the LHC and a future linear e+e^ col- 
lider. 

In the second contribution, a study of the measure- 
ment of the forward-backward asymmetry, ApB, in 
e^e~ and fi'^fJ-~ events is presented. The forward- 



backward asymmetry of events in R\m II can 

yield a measurement of the effective weak mixing angle 
sin^ ^vv and can provide a test of the standard model 
7*/Z interference at £'^£~ invariant masses well above 
the 200 GeV center of mass energy of the LEP col- 
lider. The asymmetry at large partonic center of mass 
energies can also be used to study the properties of 
possible new neutral gauge bosons, and to search for 
compositeness and large extra dimensions. Estimates 
of the statistical and systematic uncertainties expected 
in Run II for ApB and sin^ are given. The un- 
certainty for sin^ 6w is compared with the precision 
expected from LHC experiments, and from a linear 
collider operating at the Z pole. 

The third contribution summarizes the features of 
the FORTRAN package GAPP which performs a fit to 
the electroweak observables and extracts the probabil- 
ity distribution function of Mjj. 
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We discuss the prospects for measuring the W mass and width in Run II. The basic techniques used to measure Mw are 
described and the statistical, theoretical and detector-related uncertainties are discussed in detail. Alternative methods of 
measuring the W mass at the Tevatron and the prospects for Mw measurements at other colliders are also described. 



1. Introduction 

Measuring the W mass, Mw, and width, Tw are 
important objectives for the Tevatron experiments in 
Run II. The goal for the W mass measurement at the 
Tevatron in Run II is determined by three factors: 
the direct measurement of the LEP II experiments, 
the indirect determination from within the Standard 
Model (SM), and the ultimate precision on the mea- 
sured top quark mass. The expectations for LEP II 
appear to be an overall uncertainty of approximately 
±35 MeV/c2 [ 0]. The indirect determination is at 
the ±30 MeV/c^ level and is not likely to significantly 
improve given the end of the LEP and SLC programs. 
Finally, the top quark mass precision may reach the 
±2 GeV/c^ level, which corresponds to a parametric 
uncertainty of SMw = 12 MeV/c^ [||]. The constraint 
provided by a successful ± 30 MeV/c^ W boson mass 
measurement per experiment per channel]^ would have 
an impact on electroweak global fitting comparable 
to that of the LEP Z asymmetries. If the ultimate 
precision on the M\y determination could reach ~ 
±30 MeV/c^, then the bound on the Higgs boson mass 
would reach approximately ±30 GeV/c^ [ |^. With 
the best fit central value close to the current LEP II 
direct search lower hmit of Mh > 113.2 GeV/c^ [|], 



considerable pressure can be brought to bear on the 
SM. 

This document is structured as follows. The basic 
techniques used to measure the W mass and width 
are briefly reviewed in section |[ The statistical 
and detector-related uncertainties affecting the W 
mass and width measurements are discussed in more 
detail in section |^ and section 0, respectively. A 
number of systematic uncertainties clearly do not scale 
statistically and these are addressed separately in 
section |[ The expected errors on the measured W 
mass in Run II using the conventional transverse mass 
method and the W width are summarized in section ^ 
Alternative methods of measuring the W mass at the 
Tevatron are described in section |^ and prospects for 
M\Y measurements at other colliders are discussed in 
section 0. Finally, some theoretical considerations im- 
portant for future W mass measurements are brought 
up in section 0. Section 10 concludes this document. 



^ While the measurements from the different channels and 
different experiments provide cross checks, the combined mea- 
surement is not expected to yield a much better precision than 
a single measurement because of large common uncertainties. 



2. Mw and Tw Measurements from the Mt 
Lineshape 

The determination of Mw depends on the two body 
nature of the W decay: W (.vg. The kinematical 
Jacobian peak and sharp edge at the value of Mw /'^ is 
easily observed in the measurement of the transverse 
momentum of either of the leptons. In practice, the 
situation is difficult due to both challenging exper- 
imental issues and the fact that phenomenological 
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Table 1 

Advantages and disadvantages to an Mw determina- 
tion via transverse mass or lepton transverse momenta. 



Measurable 


sensitivity 


resolution sensitivity 


Mt 


small 


significant 


Pt 


significant 


small 


P^T 


significant 


significant 



assumptions must be made in order to perform the 
analysis. Because the standard measurable cannot 
be written in closed form, an unbinned maximum 
likelihood calculation is required. Figure |l| shows a 
calculation of (unsmeared) with — 0; the effect 
of finite p^ ; and the inclusion of detector smearing 
effects. It is apparent that is very sensitive to the 
transverse motion of the W boson. 

Historically, precise understanding of p^ has been 
lacking, although it is currently modeled by measur- 
able parameters through the resummation formalism 
of Collins, Soper, and Sterman [ |^. For this reason, 
the transverse mass quantity was suggested [ |j and 
has been the traditional measurable. It is defined by 



Mt 



(1) 



where 4>i^,^ is the angle between the charged lepton 
and the neutrino in the transverse plane. The observ- 
ables are the lepton transverse energy or momentum 
prf and the non-lepton transverse energy u (recoil 
transverse energy against the W), from which the 
neutrino momentum pjf and the transverse mass Mt 
are derived. Figure g shows that the sensitivity of 
Mt to p^ is nearly negligible. While considerably 
more stable to the phenomenology of the production 
model, the requirement that the neutrino direction be 
accurately measured leads to a set of experimental 
requirements which are difficult in practice to control. 
So, there are different benefits and challenges among 
the direct measurements of the transverse quantities, 
pip, Pt, and Mt- Table ^ lists these relative pros and 
cons of the transverse mass and transverse momentum 
measurements. 

Both CDF and D0 have determined the W boson 
mass using the transverse mass approach. The indi- 
vidual measurements of both experiments are shown 
in Table ^ and the overall combined result is 



Mw = 80.452 ± 0.062 GeV/c^ 



(2) 



The W boson width is precisely predicted in terms 
of well- measured SM masses and coupling strengths: 



Tw — 



GfM^ 
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a,{Mw) 
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Figure 1. The effects of resolution and the finite p^ 
on p?p in W boson decay. The histogram shows 
without detector smearing and for p^ = 0. The dots 
include the effects of adding finite p'^ , while the shaded 
histogram includes the effects of detector resolutions. 
The effects are calculated for the D0 Run I detector 
resolutions. 



x(l + 0(l%)) 
2.093 ± 0.002 GeV 



(3) 



where the uncertainty is dominated by the experimen- 
tal Mw precision [ 0, || . The mass and width of the W 
boson connect both theoretically and experimentally, 
as Tw has been extracted from a lineshape analysis 
using techniques developed for the W mass measure- 
ment. Combining CDF electron and muon data from 
1994-95 yields a result with 140 MeV precision [Isl]: 



Vw = 2.04 ± 0.11 (stat) ± 0.09 (syst) GeV. 



(4) 



In this measurement, m < 20 GeV is required to 
improve the Mt resolution and to reduce backgrounds. 

Figure |^ shows the dependence of the Mt spectrum 
on Tw In the region Mt > 100 GeV/c^, the 
lineshape is sensitive to Tw but relatively insensitive 
to uncertainties in the p^ resolution. Thus, Tw is 
extracted from a fit to the region 100 GeV/c^ < Mt < 
200 GeV/c^, after signal and background templates 
are normalized to the data in the region 40 GeV/c < 
Mt < 200 GeV/c^. Figure | shows the fits to the 
CDF electron and muon data. The upper limit Mt < 
200 GeV/c^ is somewhat arbitrary. 

The measurement of Tw depends on a precise 
determination of the transverse mass lineshape. Thus, 
the same error sources contribute to both the W mass 
and width measurement. In the following we discuss 
these sources, concentrating on how they impact the 
W mass measurement. Run II projections for the 
individual uncertainties contributing to the W width 
measurement are presented in section a. 



4 



Table 2 

Tevatron results for Mw- Nw is the number of W boson events observed. Scale and non-scale systematic errors 
are listed separately. 



Experiment 


JCdt 

pb-i 


Nw 


Mw 
GeV/c2 


± stat 
GeV/c2 


± sys 
GeV/c2 


± scale 
GeV/c2 


CDF Run e 


4.4 


1130 


79.91 


0.35 


0.24 


0.19 


CDF Run 


4.4 


592 


79.90 


0.53 


0.32 


0.08 


CDF Run la e 


18.2 


5718 


80.490 


0.145 


0.130 


0.120 


D0 Run la e 


12.8 


5982 


80.350 


0.140 


0.165 


0.160 


CDF Rim la fi 


19.7 


3268 


80.310 


0.205 


0.120 


0.050 


CDF Run lb e 


84 


30,100 


80.473 


0.065 


0.054 


0.075 


D0 Run lb e 


82 


28,323 


80.440 


0.070 


0.070 


0.065 


D0 Run lb e, forward 


82 


11,089 


80.757 


0.107 


0.091 


0.181 


CDF Run lb 


80 


14,700 


80.465 


0.100 


0.057 


0.085 




55 60 65 70 75 80 85 90 95 

(GeV) 



Figure 2. The effects of resolution and the finite 
Pjf on AIt in W ^ ev. The histogram shows Mt 
without detector smearing and for — 0. The dots 
include the effects of adding finite p^, while the shaded 
histogram includes the effects of detector resolutions. 
The effects are calculated for the D0 Run I detector 
resolutions. 
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3. Statistical Uncertainties in the Mw Determi- 
nation 

In order to reach the target precision for Mw, con- 
siderable luminosity will be required. Presuming that 
Run II is to deliver an integrated luminosity of 2 fb""'^, 
the statistical precision on Mw can be estimated from 
the existing data (see Table ||). Figure || shows the 
W statistical uncertainties in these measurements as 
a function of \/^/Nw, demonstrating a predictable 
extrapolation to Nw — 700, 000 which corresponds 
to a Run II dataset per experiment per channel. 
The statistical uncertainty from this extrapolation is 
approximately 13 MeV/c^. For a goal of ±30 MeV/c^ 
overall uncertainty, this leaves 27 MeV/c^ available in 
the error budget which must be accounted for by all 
systematic uncertainties. 



Figure 3. Dependence of the Mt spectrum on Tw- 
Each spectrum is normalized to unit area. 



4. Detector-specific Uncertainties in the Mw 
Determination 

After the lepton energy and momentum scales, the 
modeling of the W recoil provided the largest system- 
atic uncertainty in the CDF Run lb W mass measure- 
ment. Since Z statistics dominates this number, it 
can be expected to be reduced significantly in Run II. 
Non-Z related recoil systematics were estimated to 
enter at the 10 MeV/c^ level, which is probably 
indicative of the limiting size of this error. The increase 
in the average number of overlapping minimum bias 
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statistical Uncertainties 




0.02 



1/sqrt (#W) 



Figure 4. CDF 1994-95 e and /i data, on a semiloga- 
rithmic scale, with best fits for Tw- The background 
estimates are also shown. The insets show the Jacobian 
peak regions on a linear scale. 



Figure 5. Statistical uncertainties in Run 1 Mw 
measurements. Each circle represents either a CDF 
or D0 measurement. The result of a straight line 
fit is shown. The shaded box is the approximate 
extrapolation to a 2 fb~^ Run II result. 



events in Run II may seriously impact the recoil model 
systematics, although various detector improvements 
may partly compensate for this. 

Much of the understanding of experimental system- 
atics comes from a detailed study of the Z bosons and 
hence as luminosity improves, systematic uncertainties 
should diminish in kind. Certainly, the scale and 
resolution of the recoil energy against the W come from 
measurements of the Z system. Likewise, background 
determination, underlying event studies, and selection 
biases depend critically, but not exclusively, on Z 
boson data. Most importantly, the lepton energy 
and momentum scales depend solely on the Z boson 
datasets. 

Figure ^ shows the CDF and D0 systematic uncer- 
tainties for both electrons and muons as a function 
of 1/ \/Nw and in particular the calorimeter scale 
uncertainties for electrons. This latter important 
energy scale determination is currently tied to the 
determination of fiducial di- lepton decay resonances, 
notably the Z boson, but also the J/ip, T and the 
E/p dependence on the energy E, using electrons 
from W and Z decays. As the statistical precision 
improves, the dominant feature of the scale deter- 
mination becomes its value in the region of Mw, 



so offsets and any low energy nonlinearities become 
relatively less important and hence reliance on the low 
mass resonances is reduced. On the other hand, for 
the muon momentum scale determination, where the 
observable is the curvature, low mass resonances are 
also important. Figure || suggests that this uncer- 
tainty is truly statistical in nature and extrapolates 
to approximately the 15 MeV/c^ level. The ability to 
bound non-linearities using collider data may become a 
limiting source of error in Run II. Hence, the remaining 
systematic uncertainties must be controlled to a level of 
approximately 22 MeV/c^ in order to reach the overall 
goal of ±30 MeV/c^. 

Figure ^ also shows the non-scale systematic uncer- 
tainties from both the CDF and D0 electron mea- 
surements of Mw and the CDF muon measurement. 
Here the extrapolation is not as straightforward, but 
there is clearly a distinct statistical nature to these 
errors. That they appear to extrapolate to negative 
values suggests that the systematic uncertainties may 
contain a statistically independent component for both 
the muon and the electron analyses. 

For both Mw and Tw analyses, the Z ^ ££ data 
constrain both the lepton scales and resolutions and 
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Figure 6. Systematic uncertainties for each Run I 
Mm measurement. The open squares are the four 
electron measurements from CDF and D0, the circles 
are the scale uncertainties from two D0 electron 
measurements and the Run lb CDF measurement, 
and the diamonds are the systematic uncertainties 
(excluding scale) for the CDF muon measurements. 
The large box is the position of the extrapolated 
statistical uncertainties to the Run II luminosity. The 
lines are linear fits to each set of points. 



an empirical model of the hadronic recoil measurement. 
QED corrections are an issue in measuring the Z mass, 
and the discussion of these corrections should be in 
terms of the WjZ mass ratio. In a high-precision 
width measurement, more effort will also be needed 
to place bounds on possible tails in the Icpton and 
recoil resolution functions. Uncertainty in the recoil 
measurement is predominantly statistical in how well 
model parameters are determined. Several cross checks 
which improve with statistics independently ensure the 
efficacy of the model. 

Selection biases can be studied with various control 
samples, notably the second lepton originating from 
Z decays. The QCD background can be also studied 
by varying cuts and studying control samples. The 
background from W ^ tv va well understood, and the 
background from Z ^ U will be reduced for Run II 
since the tracking and muon coverages are improved 
for both experiments. 



5. Theoretical Uncertainties in the M-^ Deter- 
mination 

The A/t lineshape simulation requires a theoretical 
model, as a function of F^y and Mw, of ^. dy'dpT ' 
including correlations between pT and s. For pro- 
ducing high-statistics fitting templates, a weighted 
Monte Carlo generator is useful, so that Mw, ^w, and 
the PT spectrum can be varied simply by reweighting 
events. Because the measurement of the recoil energy 
against the W , u, is modeled empirically, the generator 
does not have to describe the recoil energy at the 
particle level. A detailed description of final-state 
QED radiation is important, because bremsstrahlung 
affects the isolation variables needed to select a clean 
W sample. 

The W and Z pr spectra are not calculable using 
perturbation theory at low px- In this region, the 
perturbative calculation must be augmented by a 
non-perturbative contribution which depends on three 



parameters (see section 5.2.1) which are tuned to fit 
the Z ^ ££ data. Theoretical guidance is useful for 
choosing an appropriate set of parameters to vary. A 
strategy such as has been used in the CDF Run lb 
analysis to use theory to extrapolate from the Z pT 
distribution to the W pt distribution seems to limit 
the effect of theoretical assumptions to ±5 MeV/c^. 

The parameters of parton distribution functions are 
also empirical, and seldom have quoted uncertainties. 
PDF uncertainties seem under control for Run I data 
but will need improvement to avoid becoming dom- 
inant in Run II. More work is needed to determine 
how both to minimize the impact of PDF uncertainties 
(e.g. by extending the lepton rapidity coverage of the 
measurements as done in the D0 analysis [ |l^ ) and to 
evaluate the effects of PDF uncertainties in precision 
measurements . 

To date, ad hoc event generators have been used in 
the W mass and width measurements. In Run II, these 
measurements will reach a precision of tens of MeV/c^, 
requiring much more attention to detail in Monte Carlo 
calculations. Precision electroweak measurements in 
Run II should strive to use (possibly to develop) pub- 
lished, well documented Monte Carlo programs that 
are common to both collider experiments. In particu- 
lar, the Mw and Tw measurements would benefit from 
a unified generator that incorporates state-of-the-art 
QED and electroweak calculations, uses a boson pT 
model tunable to Run II data, and correctly handles 
W bosons that are produced far off-shell. 

The W width uncertainty in the Mw measurement 
could become significant but assuming the SM Mw- 
Tw relation, it won't be a dominant uncertainty. 
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5.1. Parton Distribution Functions 

The transverse mass distribution is invariant under 
the longitudinal boost of the W boson. However, the 
incomplete i] coverage of the detectors introduces a 
dependence of the measured Mt distribution on the 
longitudinal momentum distribution of the produced 
VF's, determined by the PDF's. Therefore, quantifying 
the uncertainties in PDFs and the resulting uncertain- 
ties in the W mass measurement is crucial. 

5.1.1. Constraining PDFs from the Tevatron 
data 

The measurement of the W charge asymmetry at 
the Tevatron, which is sensitive to the ratio of d to 
u quark densities in the proton, is of direct benefit in 
constraining PDF effects in the W mass measurement. 
This has been demonstrated by the CDF experiment. 
Following Ref. 
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they made parametric modifica- 
tions to the MRS family of PDFs. These modifications 
with retuned parameters are listed in Table |^ and 
their predictions are compared to the W lepton charge 
asymmetry measurement and the NMC d/u data [ 
in Fig. 1^. From the variation among the six reference 
PDFs, an uncertainty of 15 MeV/c^ was taken which 
is common to the electron and muon analyses. 

Since the Run lb charge asymmetry data is domi- 
nated by statistical uncertainties, we expect a smaller 
uncertainty for the Run II measurement. Measure- 
ments of Drell-Yan production at the Tevatron can be 
used to get further constraints on PDFs. 



Table 3 

Reference PDFs and modifications 



PDF set 


Modification 


MRS-T 
MRS-R2 

MRS-Rl 


d/u d/u X (1.07 - 0.07e-**"=) 
d/u d/u + O.llx X (1 -1- x) 
d/u ^ d/u X {1.00- OMe 2(* 0.015*)') 




Figure 7. (a) The CDF measurement of the W lepton 
charge asymmetry, (b) The NMC d/u data evolved 
to = M^. The gray bands represent the range 
spanned by the six reference PDFs considered. 



5.1.3. A Global Approach 

There has been a systematic effort to map out 
the uncertainties allowed by available experimental 
constraints, both on the PDFs themselves and on 
physical observables derived from them. This approach 
will provide a more reliable estimate and may be the 
best course of action for precision measurements such 
as the W mass or the W production cross section. This 
has been emphasized at this workshop by the Parton 
Distributions Working Group [ 0| . 



5.1.2. Reducing the PDF uncertainty with a 
larger t] coverage 

Since the PDF uncertainty comes from the finite r/ 
coverage of the detectors, it is expected to decrease 
with the more complete rapidity coverage of the Run II 
detectors. The advantage of a larger rapidity coverage 
has been demonstrated by the D0 experiment: the 
uncertainty on the W mass measurement using their 
central calorimeter was 11 MeV/c^, while that using 
both the central and end calorimeter was 7 MeV/c^. 
With the upgraded calorimeters and trackers for the 
range \ri\ > 1, the CDF experiment can measure the 
W mass over a larger rapidity range in Run II. 



5.2. W Boson Transverse Momentum 

The neutrino transverse momentum is estimated by 
combining the measured lepton transverse momentum 
and the W recoil: pJf = —{p^ + u). It is clear 
therefore that an understanding of both the underlying 
W boson transverse momentum distribution and the 
corresponding detector response, usually called the 
recoil model, is crucial for a precision W mass measure- 
ment. For the CDF Run lb W mass measurement, the 
systematic uncertainties from these two sources were 
estimated to be 15 - 20 MeV/c^ and 35 - 40 MeV/c^, 
respectively, in each channel [ |l^ . 
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5.2.1. Extracting the pj^ Distribution 

The strategy employed in Run I, which is expected 
to be used also in Run II, is to extract the underlying 
distribution from the measured p^ distribution {Y 
is the weak boson rapidity): 



(fa (fa/dp)^dY 



dp^dV dp^dY d'^a/dp^dY 



(5) 



where the ratio of the W and Z differential distri- 
butions is obtained from theory. This method relies 
on the fact that the observed distribution suffers 
relatively little from detector smearing effects, allowing 
fits to be performed for the true distribution. The 
CDF Run lb data and the results of a Monte Carlo 
simulation using the best fit parameters are compared 
in Fig. I 
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Figure 8. The observed distributions in the 
electron and muon channels using the CDF Run lb 
data. Also shown are the curves for the Monte Carlo 
simulation using the best fit parameters for the input 
distribution. 



The experimental uncertainties, as in many as- 
pects of the W mass measurement, are dominated 
by the available Z statistics and should 

scale correspondingly with the delivered luminosity in 
Run II. Theoretical uncertainties in the ratio of W 



to Z transverse momentum distributions contribute 
a further 0(5) MeV/c^ to the overall error. The 
two sources of uncertainty are compared for the CDF 
Run Ih W fiiy analysis in Fig. ^. 
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Figure 9. (a) A comparison of the two sources of 
uncertainty on the derived p^ distribution, (b) The 
Pj^ distribution extracted for the CDF Run lb ^ 
fiv mass measurement. 



The ratio of W to Z transverse momentum distri- 
butions used in Eq. (^) is taken from resummation 
calculations, which attempt to resum terms corre- 
sponding to multiple soft and coUinear gluon emission 
to all orders. They thereby include the dominant 
contribution to the cross section at small boson pT 
that is missing in fixed order calculations. These 
perturbative calculations need to be augmented by a 
non-perturbative contribution which, in the case of 
impact-parameter (6) space resummations, is typically 
parameterized as a Sudakov form-factor with the fol- 
lowing form: 



F 



- exp [ -g,b^ - g^h" HQ/2Qo) 
-gig3&ln(100xiX2)] , 



(6) 



where Qo is a low scale of O(few) GeV and the 
parameters gi , 52 and g^ must be obtained from fits 
to the data [ |l5|. D0 has shown that the Run I p^ 
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data is as sensitive to gi and 32 as the low-energy 
Drell-Yan data that has largely been used to constrain 
these parameters in the past [ |l6|. The Run II data 
will therefore provide significant new constraints on the 
form of the non-perturbative contribution to the pj^ 
distribution. 

Moreover, recent theoretical developments in com- 
bining the advantages of 6— space and pt— space resum- 
mation formalisms may provide a better theoretical 
framework for extracting the underlying distribu- 
tion in Run II [ |l|] . 

In short, the precision Z data available in Run II 
together with further theoretical advances will reduce 
in a number of ways the systematic uncertainties due 
to the knowledge of the p}^ distribution, perhaps down 
to the level of 5 MeV/c^. 

5.3. QCD Higher Order Effects 

The W bosons are treated as spin-one particles and 
decay via the weak interaction into a charged lepton 
(e, /i or t) and a neutrino. The charged leptons 
are produced from the W decay with an angular 
distribution determined by the Olaj.) calculation of [ 
which, for W'^ bosons with a helicity of -1 with 
respect to the proton direction, has the form : 

-j-^ — oc 1 -f aiipr) cos 6*03 + 02 (pt) cos^ Ocs (7) 
a cos Ocs 

where px is the transverse momentum of the W and 
6*03 is the polar direction of the charged lepton with 
respect to the proton direction in the CoUins-Soper 
frame [ |9|. ai and 02 are px dependent parameters. 
For pt — 0, ai — 2 and 02 = 1, providing the 
angular distribution of a boson fully polarized along 
the proton direction. For the p}^ values relevant to 
the W mass analysis {p^ < 30 GeV/c), the change 
in W polarization as increases only causes a 
modest change in the angular distribution of the decay 
leptons 

While the uncertainty associated with the change in 
the angular distribution of the W decay lepton due 
to higher order QCD corrections (a few MeV/c^) has 
been negligible for the Run I measurements, it can 
not be ignored for the Run II measurements (see the 
Photon and Weak Boson Physics working group report 
for more details). 

5.4. QED Radiative Effects 
5.4.1. Introduction 

The understanding of QED radiative corrections is 
crucial for precision W mass measurements at the 
Tevatron. The dominant process is final state radiation 
(FSR) from the charged lepton, the effect of which 
strongly depends on the lepton identification criteria 
and the energy or momentum measurement methods 



employed. Calorimetric energy measurements, such 
as those employed in the electron channel, are more 
inclusive than track based momentum measurements 
used in the muon channel and the effect of FSR 
is consequently reduced. In the CDF Run lb W 
mass measurement the mass shifts due to radiative 
effects were estimated to be —65 ± 20 MeV/c^ and 
-168± 10 MeV/c^ for the electron and muon channels, 
respectively [ These effects will be larger in 

Run II due to increase in tracker material in CDF and 
magnetic tracking in D0. 

The Monte Carlo program used for the Run I W 
mass measurement incorporated a calculation of QED 
corrections by Berends and Kleiss [ This treat- 
ment, however, does not include initial state radiation 
(ISR) and has a maximum of one final state photon. 
The effect of multiple photon emission was estimated 
by comparing the calculation of Berends and Kleiss 
to PHOTOS [ |l|, a universal Monte Carlo program 
for QED radiative corrections that can generate a 
maximum of two final state photons. Likewise, the 
effect of ISR and other missing diagrams was estimated 
by comparing the calculation of Berends and Kleiss 
to a full 0{a) calculation by Baur et at [ |2^. The 
resulting systematic uncertainties on the W mass are 
estimated to be 20 MeV/c^ and 10 MeV/c^ in the 
electron and muon channels, respectively [|l4j. Clearly 
these systematic uncertainties become much more 
significant in the context of statistical uncertainties of 
0(10) McV/c^ expected for 2 fh'^ in Run II. 

The next section describes in more detail the calcu- 
lation by Baur et ai, which forms the basis for a new 
event generator. Some studies of the effects of QED 
radiation on the W mass measu rement are presented 
in section 5.4.4 . Section 5.4.5 briefiy outlines some 
work in progress that should further reduce systematic 
uncertainties due to radiative corrections in Run II. 



5.4.2. WGRAD 

WGRAD is a program for calculating 0{a) electroweak 
radiative corrections to the process qq' ^^z/, 
including the real photon contribution qq' — > 
Both ISR from the incoming quarks, FSR from the 
final state charged lepton, and interference terms are 
included. Many more details can be found in [ ^ . 

The most important generator level cuts are on 
the final state photon energy and coUinearity for 
radiative events. The photon energy cut, controlled 
by the parameter Ss, is made on the fraction of 
the parton's energy carried by the emitted photon 
in the parton-parton center of mass system: E* > 

6s^/S/2. The photon coUinearity cut, controlled by 
the parameter Sc, is made on the angle between the 
charged fermion and the emitted photon in the same 
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Table 4 

The fraction of qq' — s- iJ-i^i"/) events containing a final 
state photon for different final state photon soft and 
coUinear cuts. Events are generated with ISR only, 
FSR only, and with a full treatment of QED radiation. 



Photon Cuts 


ISR 


FSR 


Full 


Photon Cuts 


ISR 


FSR 


Full 


Ss = 0.01, (5c = 0.01 


1.6% 


9.4% 


11.1% 


Ss = 0.01, (5c = 0.001 


2.5% 


9.4% 


12.0% 


(5, = 0.001,5c = 0.001 


4.1% 


15.5% 


20.0% 


6s = 0.001, (5c = 0.0001 


5.2% 


15.5% 


21.3% 



frame: cosd* < 1 — (5c. However, final state coUinear 
singularities are regulated by the finite lepton masses 
and the above cut is only implemented for quarkonic 
radiation when ISR is included. The fraction of 
radiative events corresponding to different photon cuts 
is given in Table ^ for the process qq' l^i^il) at 
\/s = 2 TeV. Loose fiducial cuts p!^ > 10 GeV/c, 
|77^| < 2 and > 10 GeV/c have been apphed. The 
inclusion of ISR increases the photon yield by around 
30%, depending on the soft and coUinear photon cuts 
applied. The fractions are significantly higher for the 
process qq' in the cases that FSR is included. 

The effect on the fitted W mass of the inclusion of ISR 



is examined in section 5.4.4 



5.4.3. Event Generation 

WGRAD has been turned into an event generator 
through a suitable unweighting scheme described ex- 
tensively in [ |2^. A significant complication is the 
presence of negative qq' t^y event weights in 
the program which, while expected to cancel with 
positive qq' — > i^v^f event weights in the calculation 
of physical observables, nevertheless appear separately 
in the unweighting procedure. The approach here 
has been to unweight the negative weight events in a 
similar manner to the positive weight events, such that 
the output consists of both positive and negative unit 
weight events. The fraction of negative weight events, 
plotted in Fig. ^for the process qq' — > ^J-iy{"f), depends 
strongly on the soft and coUinear photon cuts applied. 
It is not significantly different for qq' — > ev^j) events. 
The effect of negative weights on the fitted W mass is 
examined in the next section. 



83 = 0.001, 5^ = 0.0001 
83 = 0.001, 8^ = 0.001 
8 =0.01, 8 =0.001 
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Figure 10. The negative weight fraction versus for 
different soft and coUinear photon cuts. 



5.4.4. The Effect of QED Radiation on the 
Measured W Mass 

WGRAD has been used to generate large W ^ ^iv 
event samples for the purposes of investigating the 
effect of QED radiation on the measurement of the 
W mass. The events have been generated at = 
1.8 TeV in order to make use of the CDF Run I 
W production model and detector smearing param- 
eterizations. The W production model, extracted 
from the Run I Drell-Yan data, is used to smear 
the true W transverse momentum. The CDF recoil 
model is then used to translate this into a measured 
Pt{W), which is combined with the smeared lepton 
and photon momenta to form a realistic transverse 
mass distribution. Loose fiducial cuts > 20 GeV/c, 
Iry'^l < 2 and missing- _Et > 20 GeV are apphed. To 
simulate the CDF muon identification criteria, events 
are rejected if a photon with > 2.0 GeV is found 
within an r] — (j) cone of radius 0.25 around the muon. 
Low energy photons inside the cone are not included 
in the measurement of the muon pT, as is the case 
experimentally. 

The unweighted event samples, all generated with 
Mw = 80.4 GeV/c2 and F^ = 2.1 GeV, are divided 
into "data" and "Monte Carlo" sub-samples and fitted 
against one another in pseudo-experiments. The fit is 
to the transverse mass distribution in the range 50 < 
Mt < 100 GeV/c^. For a number of events in the 
transverse mass fit region equal to that in the CDF 
Run IhW-^fiv analysis, the resulting statistical error 
is very similar. 

As a cross check of this procedure, "data" and 
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Table 5 

The results of pseudo-experiments generated by WGRAD with different treatments of QED radiative effects. Further 
details are given in the text. 





"Data" 


"Monte Carlo" 


Fit Result 


(a) 


FULL; Ss = 0.01, (5c = 0.01 


FULL; Ss = 0.01, (5c = 0.01 


80.4030 ±0.0069 GeV/c^ 


(b) 


FULL; = 0.01,5c = 0.01 


FULL no neg.; Ss = 0.01, 6c = 0.01 


80.4027 ±0.0069 GeV/c^ 


(c) 


FSR only; 6^ = 0.01, (5c = 0.01 


FULL; Ss = 0.01, Sc = 0.01 


80.392 ±0.006 GeV/c^ 


(d) 


FSR only; (5, = 0.001, 6c = 0.001 


FULL; Ss = 0.01, (5c = 0.01 


80.381 ±0.006 GeV/c2 


(e) 


FULL; (5, = 0.001, (5c = 0.001 


FULL; ,5, = 0.01, (5c = 0.01 


80.389 ±0.009 GeV/c^ 


(f) 


FULL no neg.; Ss = 0.001, S^ = 0.001 


FULL; 6s = 0.01, ,5c = 0.001 


80.399 ±0.009 GeV/c^ 



"Monte Carlo" samples generated with identical cuts 
are fitted against one another, with the result shown 
in Table ||(a). It is interesting to note that if the 
negative weight events, which occur at the 0.2% level in 
the "Monte Carlo" sample, are removed, the fit result 
changes by less than 0.5 MeV/c^ (Table |(b)). 

Table ||(c) shows the result of fitting "data" gener- 
ated with FSR only. The shift in the fitted W mass 
of « 8 MeV/c^ is consistent with the estimate given 
in [ Q of the effect of ISR on the fitted W mass, 
although the uncertainties here are rather large. If 
the soft and collinear cuts are reduced in the "data" 
sample, as shown in Table ^(d), the fitted W mass 
shifts significantly downwards. This is to be expected 
since the track based muon px measurement does 
not incorporate collinear photons. The setting of 
soft and collinear photon cuts is therefore particularly 
important in the generation oi W ^ Monte Carlo 
samples. 

The fits shown in Table ||(e) and (f) are performed 
in order to examine the effect of negative weights on 
the fit when, as in the case of this "data" sample, 
negative weights are present at the 5% level. When 
the negative weight events are excluded from the fit, 
the result changes by 10 MeV/c^. The larger shift 
in the fitted W mass with respect to Table ||(b) is 
commensurate with the larger negative weight fraction 
in this sample. 

5.4.5. Work in Progress 

A remaining source of systematic uncertainty due 
to QED radiation is the effect of multiple photon 
emission. As discussed above, this has previously been 
estimated by comparing the Berends and Kleiss single 
photon calculation with the results of running the 
PHOTOS algorithm. Recently, however, complete matrix 
element calculations of the processes qq' £='=1/77 and 



qq £+£^77 have been performed [ It may be 
possible in the future to do detailed comparisons of the 
results of these calculations and the PHOTOS algorithm, 
in order to arrive at a better constrained systematic 
uncertainty due to multiple photon emission. 

Furthermore, a complete set of 0{a) electroweak 
radiative corrections to the process qq £^£^, 
including the real photon contribution qq £+£^7, 
will soon be available. This will enable a consistent 
Monte Carlo description of the W data and the Z data, 
upon which the W mass analysis crucially depends for 
the understanding of gauge boson production and the 
calibration of the detectors. 

5.4.6. Summary and Conclusions 

Systematic uncertainties due to QED radiative ef- 
fects currently run at the level of w 20 MeV/c^ in 
the electron channel and « 10 MeV/c^ in the muon 
channel. A large contribution to this uncertainty is 
the effect of ISR and interference terms, which are not 
present in the Berends and Kleiss calculation and the 
PHOTOS algorithm that have previously been used in W 
production Monte Carlo programs. 

A full 0{a) calculation by Baur et al. has been 
used as the basis for a new event generator. The 
results of several pseudo-experiments generated with 
different treatments of QED radiative effects agree 
with previous estimates. They show that negative 
weights need to be treated carefully, especially in the 
case of very small soft and collinear photon cuts. 

Further studies of QED radiative corrections to W 
production will continue as new calculations become 
available. It is clear, however, that the use of new 
programs such as WGRAD could significantly reduce 
systematic uncertainties due to QED radiative cor- 
rections in Run II, either through explicit corrections 
being applied to the extracted W mass, or through 



12 



their use in new Monte Carlo event generators. The 
remaining systematic uncertainties due to QED correc- 
tions might then be reduced to the level of 5 MeV/c^ 
and 10 MeV/c^ in the muon and electron channels, 
respectively. 

6. Summary of Run II Expectations 

As has been discussed in previous sections, many 
of the systematic uncertainties in the W mass mea- 
surement approximately scale with statistics. These 
are listed in Table ^ for the Run lb CDF muon 
analysis and should scale to w 20 MeV/c^ for an 
integrated luminosity of 2 fb^^. With reasonable 



Table 8 

Dominant uncertainties for contrasting components of 
the D0 Mw determination. The quantities shown 
are the shift in Mw for a la change in the relevent 
parameter. The EM resolution term refers to the 
sampling term for the resolution function. Taken from 
Ref. [|2l. 



Source 


SMw{Mt) 


5Mw{p't) 


w 
Pt 


10 


50 


EM resolution 


23 


14 


hadron scale 


20 


16 


hadron resolution 


25 


10 


backgrounds 


10 


20 



Table 6 

Errors on the CDF Run lb muon W mass which scale 
statistically, in MeV/c^. 



Source 


error 


Fit statistics 


100 


Recoil model 


35 


Momentum resolution 


20 


Selection bias 


18 


Background 


25 


Momentum scale 


85 



assumptions for the size of non-scaling systematics 
such as those due to PDFs and higher order QED 
effects, a 40 MeV/c^ measurement in the muon channel 
by each experiment seems achievable. The systematic 
uncertainties in the electron channel are less easy to 
extrapolate given the particular sensitivity to calorime- 
ter scale non-uniformities in this channel and the extra 
material in the Run II tracking detectors. The detailed 
understanding of detector performance is of course 
difficult to anticipate, although it is clear that both 
scalable and non-scaling systematics would be easier 
to understand if fast Monte Carlo generators including 
all the relevant effects were available. 

The individual uncertainties for the Run lb Ty/ 
measurement are listed in Table ^ together with their 
projections for 2 fb~^. All but the last three sources 
of error are constrained directly from collider data, 
and hence should scale roug hly as 1/VZ. While the 
last three uncertainties may decrease somewhat as new 
measurements and calculations become available, they 
will not scale statistically with the Run II dataset. 
Assuming no improvement in these three uncertainties, 
while all others scale statistically, each experiment can 



make a ^40 MeV width measurement, combining e 
and ^ channels for a 2 fb~^ dataset. 

7. Other Methods of Determining Mw at the 
Tevatron 

While the traditional transverse mass determination 
has been the optimal technique for the extraction of 
Mw in the low-luminosity running at hadron colliders, 
other techniques have been or may be employed in the 
future. These methods may shuffle or cancel some of 
the systematic and statistical uncertainties resulting in 
more precise measurements. 

7.1. Transverse Momentum Fitting 

As noted above, the most obvious extensions of the 
traditional transverse mass approach to determining 
Mw are fits of the Jacobian kinematical edge from 
the transverse momentum of both Icptons. D0 has 
measured Mw using all three distributions and the 
uncertainties are indeed ordered as one would expect: 
The fractional uncertainties on Mw from the D0 
Run I measurements for the three methods of fitting 
are: 0.12% (Mt), 0.15% (pf,), and 0.21% (p^). As 
expected, the method is slightly less precise than 
the transverse mass. However, for a central electron 
(|?7| < 1), the uncertainty in the pf^ measurement 
due to the model is 5 times that in the Mt 
measurement. As can be seen from Table ||, this is 
nearly balanced by effects from electron and hadron 
response and resolutions which are relatively worse for 
Mt- Accordingly, when there are sufficient statistics 
to enable cuts on the measured hadronic recoil, the 
measurement uncertainty from the model might be 
better controlled and enable the measurement to 
compete favorably with the Mt measurement which 
relies so heavily on modeling of the hadronic recoil. 
In order to optimize the advantages of all three mea- 
surements, the D0 final Run I determination of M-w 
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Table 7 

Sources of error for CDF 1994-95 Tw measurement and extrapolations to 2 fb^^. The last three uncertainties are 
common to the e and fi analyses. 





CDF 1994-95 (-^ 2fb"') 


Source 


AF (e, MeV) 


AF in, MeV) 


Statistics 


125 (^ 30) 


195 (^ 45) 


Lepton E or px non-linearity 


60 


5 


Recoil model 


60 


90 


W Pt 


55 


70 


Backgrounds 


30 


50 


Detector modeling, lepton ID 


30 


40 


Lepton E or px scale 


20 


15 


Lepton resolution 


10 


20 


PDFs (common) 


15 


15 


Mw (common) 


10 


10 


QED (common) 


10 


10 


Uncorrelated systematic 


112 (^ 25) 


133 (^ 30) 


Correlated systematic 


21 


21 


Total systematic 


115 (^ 33) 


135 (^ 37) 


Total stat -1- syst 


170 (^ 45) 


235 (^ 60) 



combined the separate results [ 25 . 

The resolution sensitivity for muon measurements is 
even less than that for electrons so that has the benefit 
of slightly favoring a transverse mass measurement 
with muons over that for electrons. 

7.2. Ratio Method 

D0 has preliminarily determined Mw by consideration 
of ratios of W and Z boson distributions which are 
correlated with Mw [ The principle is that one 
can cancel common scale factors in ratios and directly 
determine the quantity r""^"* = which can be 

compared with the precise LEP M§^^ . The quantities 
that have been considered are: 

1. r{MT) and r{pT), which has the advantage of 
being well-studied [ There are challenges 
with this approach which will be discussed below. 

2. r{E'^) which has the advantage that the peak of 
the distribution is precisely correlated with Mw, 
but the disadvantage that statistical uncertainty 
washes out the position of that peak. 

3. The difference of transverse mass distributions 
(not as precise as ratios). 

The procedure is to compare two distributions, one for 
W bosons and a similarly constructed one for Z bosons, 
for example, f^'^{x) as a function of a given variable, 
such as a; = Mt ot x = p^. Practically speaking, the 
Z boson decay electrons are scaled by a factor s and 
f^{x,s) is compared with f^{x) as a function of x, 
for different trial values of s. A statistical measure 



(the Kolmogorov-Smirnov test) is calculated for each 
s and the value of the highest Kolmogorov-Smirnov 
probability, s*"^^*, is declared to be r""^"^ and the 
desired mass is then extracted from Mw = r™*^"* x 
M^^^ . In principle, minimal Monte Carlo fitting is 
required, as the measurement is performed with data. 

Figure ^ shows the idea with an unsmeared Z boson 
transverse mass distribution compared to a simulated 
(unsmeared) W boson distribution. Various values of 
s lead to various mismatches between {Mt) and 
f^{Mj') which can be characterized by a Kolmogorov- 
Smirnov probability as a function of Mw = sx M^^ . 
This probability distribution for an ensemble of 100 
Monte Carlo experiments is shown in Fig. ^ resulting 
in an RMS of 40 MeV/c^. 

However, there are challenges to be faced using this 
technique. 

• Many systematic effects cancel in this method, 
such as electromagnetic scale, hadronic scale, 
angular scale, luminosity effects. However, these 
are first-order cancellations, some of which in 
the end are not sufficient: the second order 
effects from these quantities must be considered. 
Likewise, most resolutions have additive terms 
which do not cancel in a ratio. 

• The statistical precision of the Z sample is 
directly propagated into the resultant overall 
SMw, in contrast to the traditional approach 
where the Z boson statistics is a component of 
various of the measured resolutions. 

• The detector modeling must take into account 
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Figure 11. For unsmcared Monte Carlo events, the 
transverse mass of simulated W bosons (histogram) is 
overlayed with that of Z boson (dots) events in which 
the electron has been scaled by a factor which produces 
the best match. 



Figure 12. The Kolmogorov-Smirnov probability dis- 
tribution for various scales in a comparison of W 
and Z boson unsmeared events corresponding to 100 
experiments. The RMS is 40 MeV/c^ for 20000 events. 



small, but important differences between Z and 
W events such as underlying event, resolutions, 
efficiencies, acceptances, and the effects of the 
"extra" electron in Z boson events which compli- 
cates underlying event and recoil measurements. 

• From the physics model, there are also differences 
between the two samples which must be consid- 
ered, such as the fact that the production of Z 
and W bosons take place from the annihilation of 
like and unlike flavored quarks, respectively and 
that weak asymmetries lead to different decay 
angular distributions. 

• Particularly difficult is the need to "extra-smear" 
the electrons from Z boson decays. This is due 
to the fact that values for the heavier Z 
boson are harder, resulting in a different average 
resolution smearing. This same effect is true for 
the recoil distributions between Z and W bosons. 

• Finally, the acceptances for the two bosons are 
different since there are potentially two opportu- 
nities to select a Z boson event at the trigger 
and event selection stages. Similarly, there is 
an acceptance difference in the opposite direction 
due to electrons in Z bosons being lost in cracks 
between the CC and EC calorimeters in the D0 
detector. 



An analysis from Run la data from the D0 experiment 
has been done [|2^. Figure [l^ shows data for the scaled 
comparison and the unsealed original distributions. 
Electrons from the W boson events were selected to 
have p?p > 30 GeV/c, while those from Z boson 
events, must satisfy > 34.1 GeV/c. Electrons 
from the W sample and at least one electron from 
the Z samples were required to be in the central 
calorimeter. This results in 5244 W bosons and 535 Z 
boson events. Backgrounds are subtracted according 
to the traditional analysis. "Extra-smearing" is done 
for each accepted Z boson event (twice, for both 
electrons) 1000 times, using a different random seed 
for each smearing. Differences in the W and Z boson 
production mechanisms and acceptances result in an 
effective correction of 109 MeV/c^, while the difference 
in radiative corrections results in an effective correction 
of —116 MeV/c^. The magnitude of these corrections 
is not very different from corrections within the tra- 
ditional technique and the demand on knowing the 
uncertainties in them is similarly stringent. Figure ^ 
shows the probability distribution for the result. The 
preliminary result from this analysis for central. Run la 
electrons is 

Mw = 80.160 ± 0.360 ± 0.075 GeV/c^. 

Comparison with the traditional Run la result from 
the same data is readily made, but most appreciated 
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Figure 13. (a) The transverse mass distribution 
of W (solid) and scaled Z (dots) bosons is shown 
along with the hatched fit window, (b) The original 
distributions are shown, along with their difference. 
The Z distribution has been normalized to that of the 
W boson sample. 



with a slightly different interpretation of the Run la 
uncertainties. The Run la result [ |28| from Table || is 

80.350 ± 0.140 ± 0.165 ± 0.160 GeV/c^ 

where the first error is the statistical uncertainty 
(from W events), the second is the systematic un- 
certainty and the third is the electron scale deter- 
mination. It is important to note that the scale 
uncertainty is almost completely dominated by the 
Z boson statistics. Therefore, as a statistical uncer- 
tainty, it can be combined with the W uncertainty 
of 140 MeV/c^ for the purposes of comparison with 
the ratio method. This results in an overall "statis- 
tical" uncertainty of 212 MeV/c^. Now, the stronger 
systematic power of the ratio method is apparent (75 
versus 165 MeV/c^) and the poorer statistical power 
(360 versus 212 MeV/c^) is also evident. 

7.2.1. Prospects for Run II 

This apparent systematic power of the ratio method 
can only fully be realized in high luminosity running, 
such as Run II. The ratio method analysis of the 
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Figure 14. The Kolmogorov probability distribution 
(hatched) is shown as a function of the W boson 
mass used as a scale factor. The dotted curve is the 
Kolmogorov Likelihood and the dashed curve is the 
X^/ndf distribution (right axis). 



D0 Run lb data was recently completed [ |2^. The 
Run lb sample has 82 pb"^ of data (1994-1995 data 
set), 33,137 W and 4,588 Z events (electrons in 
both Central and End Calorimeters of D0) after the 
standard electron selection cuts. The W mass resulting 
from the ratio fit is Mw = 80.115 ± 0.211 (stat.) ± 
0.050 (syst.) GeV/c^. The statistical uncertainty is in 
good agreement with an ensemble study of 50 Monte 
Carlo samples of the same size (80.36 ± 0.25 GeV/c^). 

Early efforts at predicting the results for a Run II 
sample of 100,000 W bosons is shown in Fig. |l5| with 
full detector acceptances and resolutions taken into 
account. The statistical precision from this fit is of the 
order of 20 MeV/c^ and the systematic uncertainties 
may be nearly negligible. 

8. Prospects for Measuring Mw at Other Ac- 
celerators 

8.1. LEP II 

The prospects for determination of at LEP II 

have become fully understood in the last year with 
the accumulation of hundreds of pb~^ at four center 
of mass energies. Here we review the status as of 
the Winter 2000 conferences and project the prospects 
through to the completion of electron-positron running 
at CERN. For a review, see Ref. [ |o[ H. 

8.1.1. Data Accumulation 

The annihilation of e+e~ into W boson pairs occurs 
via three diagrams: a t-channel neutrino exchange and 
s-channel Z or 7 exchange. The final states from the 
decays of the two W bosons are: both W bosons decay 
into hadrons {qqqq, "4-q" mode); one W decays into 
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Figure 15. Same as Fig. |l|a, but for 100,000 W bosons. 



quarks, and the other into leptons (e, /i , t and 
their neutrinos, ^'qq£v"); and both W bosons decay 
leptonically. Collectively, the latter two modes are 
referred to as "non-q" modes. The efficiencies and 
sample purities are typically quite high, as shown in 
Table i. 



Table 9 

Efficiencies and sample purity for a representative 
LEP II experiment (OPAL). 



Channel 


e{%) 


purity (%) 


mm 


85 


80 


qqev 


85 


80 


qq^iv 


87 


80-90 


qqrv 


66 


80 



Table 10 

Approximate accumulated running per experiment. 
The 2000 totals are current as of the first week in 
July 2000. 



year 


beam energy (GeV) 


jCdt (pb-^) 


1996 


80.5-86 


25 


1997 


91-92 


75 


1998 


94.5 


200 


1999 


96-102 


250 


2000 


100-104 


100 



Table 11 

Typical systematic uncertainties on Mw for a generic 
LEP II experiment from data sets corresponding to 
the 189 GeV running. Entries are approximate and 
broad averages meant to give a relative sense of scale 
only. The Source labels are generally self-explanatory, 
with CR/BE standing for "Color Reconnection" and 
"Bose-Einstein" respectively. 



Source 


qqliy (MeV/c^) 


qqqq (MeV/c^) 


ISR 


15 


15 


frag 


25 


30 


4 fermion 


20 


20 


detector 


30 


35 


fit 


30 


30 


bias 


25 


25 


bckgrnd 


< 5 


15 


MC Stat 


10 


10 


Subtotal 


61 


67 


LEP 


17 


17 


CR/BE 




60 


Total 


63 


85 



The results by Spring 2000 come from running at 
center of mass energies: 172, 183, 189, and several 
energies between 190 GeV and 200 GeV. There is 
recent running above 200 GeV for a total of more 
than 400 pb~^ accumulated per experiment. Mw 
results are final for all four experiments for the 172 
and 183 GeV sets [ ^ [3^ and preliminary 
for the 1998 189 GeV running \m, WHM.M- In 
addition, ALEPH 
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L3 [ and OPAL [ |4J have 
preliminary results from the collection of runs in the 
range from 190 GeV to 200 GeV. Table |l^ shows the 
approximate accumulated running to date (July 2000). 

There are broadly two methods employed for de- 
termining Mw at LEP II. The first method is the 
measurement of the threshold of the WW cross section 
and the second is the set of constrained fits possible 
for the various measured final states. The latter set of 



methods constitute the prominent results and employ 
construction of invariant masses making use of the 
beam constraints. There are a variety of methods, 
some of which make use of the constraint Mwi = Mw2 
and some of which involve sophisticated multivariate 
analyses. The spirit of approach is much like the 
strategies employed in the top quark mass analyses of 
CDF and DO. 

The results are treated separately for the qqiv and 
qqqq final states due to the significant differences in 
systematic uncertainties. Typical _uncertainty contri- 
butions are listed in Table 
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J. Many of the 
experimental uncertainties, such as scale, background, 
and Monte Carlo generation, are statistically limited. 
For example, there is a fixed amount of ^/s = Mz 
running in each running period and that contributes a 
statistical component to the energy scale uncertainty. 
The dominant uncertainty comes from the final state 
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effects in the qqqq channel. Because the outgoing 
quarks can have color connections among them, the 
fragmentation of the ensemble of quarks into hadrons 
are not independent. This leads to an theoretical 
uncertainty called "Color Reconnection" (CR). In 
addition, since the hadronization regions of the W'^ 
and W~ overlap, coherence effects between identical 
low-momentum bosons originating from different Ws 
due to Bose-Einstein (BE) correlations may be present. 
The combined total of these two effects is currently 
accepted to contribute 52 MeV/c^ of uncertainty to the 
qqqq results. Ultimately, the non-CR/BE uncertainty 
will likely be the uncertainty in modeling single-quark 
fragmentation and associated QCD emission effects. 

8.1.2. Results, April 2000 

The preliminary results for AIw from the combined 
data taking through 1999 running period are shown 
in Table ^ The combined LEP result for the qq£v 
channels is [ Q : 



M, 



W 



80.398 ± 0.039 ± 0.031 ± 0.017 GeV/c^ 



where the first error is statistical, the second sys- 
tematic, and the third the LEP energy scale. The 
combined preliminary result for the qqqq channel is: 

= 80.408 ± 0.037 ± 0.031 ± 0.016 ± 0.052 GeV/c^ 

where the first three errors are the same as for the 
qqlv result and the fourth error is due to the combined 
CR/BE theoretical uncertainty. Taking into account 
the correlations, the combined preliminary result from 
constrained fitting for all channels is: 

= 80.401 ± 0.027 ± 0.031 ± 0.017 ± 0.018 GeV/c^ 

where the four errors are in the same order as for the 
qqqq result. The current overall result comes from 
combining the above with that from the threshold 
measurement of 



a(E) 
W 



80.400 ± 0.220 ± 0.025 GeV/c^ 



Here the first error is combined statistical and system- 
atic and the second error is the error due to LEP energy 
scale. This results in the preliminary overall LEP II 
(April 2000) value of 



Mw'^^ = 80.401 ± 0.048 GeV/c^ 



8.1.3. Prospects for the Future 

The current results are preliminary and running is 
underway at this writing with the end of LEP II sched- 
uled for the beginning of October, 2000. Eventually, 
the 1999 data will be fully analyzed and, with the 
accumulation of the final 2000 running, should result 



Table 12 

Preliminary LEP II results for AIw by experiment and 
according to reconstructed channel. The results are 
from the combination of 1996-1998 running (all exper- 
iments) plus preliminary results from 1999 running for 
ALEPH, DELPHI, L3, and OPAL. 





Mw (GeV/c^) 


Experiment 


qqiv 


qqqq 


ALEPH 


80.435 ±0.079 


80.467 ±0.086 


DELPHI 


80.230 ±0.140 


80.360 ±0.115 


L3 


80.282 ±0.102 


80.489 ±0.132 


OPAL 


80.483 ±0.078 


80.380 ±0.103 



in a combined statistical and systematic uncertainty 
(excluding the CR/BE and LEP contributions) of 
approximately 35 MeV/c^ [ |l|. With the overall 
contribution of 18 MeV/c^ and 17 MeV/c^ fr om the 
CR/BE and LEP errors respectively, the ultimate limit 
from LEP II W boson pair determination of Mw 
should be approximately 40 MeV/c^. 

8.2. LHC 

It was pointed out several years ago 



45 that the 



LHC has the potential to provide an even more precise 
measurement of Mw- This suggestion was based on 
the observations that the precision measurement of 
Mw at hadron colliders has been demonstrated to be 
possible; that the statistical power of the LHC dataset 
will be huge; and that triggering will not be a problem. 
These authors estimated that Mw could be determined 
to better than 15 MeV/c^. More recently, the ATLAS 
collaboration has studied the question in more detail [ 
Q and arrived at an uncertainty of 25 MeV/c^, per 
experiment, in the electron channel alone. 

Achieving such precision will require substantial 
further reduction of theoretical and systematic un- 
certainties, all of which must be reduced to the < 
10 MeV level. This includes the contributions from 
the W production model, parton distributions, and 
radiative decays, as well as experimental systematics 
such as the energy-momentum scale of the detector and 
any complications from underlying energy deposition 
even in the low luminosity running. While some 
have questioned whether such "heroic" progress will 
ever be possible, we would argue that it is futile to 
debate the question at this time. Rather, the best 
indicator of future LHC precision will be to see how 
well the Fermilab experiments manage to deal with the 
significant improvements in systematics which will be 
necessary in order to match the anticipated precision 
for mt- The point to be made is this: should it prove 
necessary to determine the W mass to a precision of 
10-20 MeV/c^, the LHC will have the statistical power 
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to continue the hadron collider measurements into this 
domain. The success of such a program will then 
depend on 

• Consensus in the field that such precision is 
needed. One such justification might be to distin- 
guish among different models of supersymmetry- 
breaking using global fits including Mw, the top 
mass and the light Higgs mass. It is likely that 
a big parallel effort to push down the top mass 
uncertainty to the 1 GeV/c^ level would also then 
be needed; 

• A major, multi-year effort within the LHC ex- 
periments in order to understand their detectors 
and their response to leptons, missing transverse 
energy and recoil hadrons at the required level. 
This is a measurement which places heavy bur- 
dens on manpower as it requires an understand- 
ing of the detector which is more precise than for 
any other measurement; 

• A comparable major effort to reduce the theo- 
retical uncertainties through better calculations, 
through control-sample measurements, and work 
on parton distributions. 

This is not a program that will be undertaken lightly. 
But should it turn out to be necessary, the experience 
of Run II at the Tevatron will be invaluable in carrying 
it out. 

8.3. A Linear Collider 

The W mass can be measured at a Linear Col- 
lider (LC) in W~^W~ production either in a dedi- 
cated threshold scan operating the machine at ^/s « 
161 GeV, or via direct reconstruction of the W bosons 
in the continuum {y/s — 0.5—1.5 TeV). Both strategies 
have been used with success at LEP II. 

In the threshold region, the W^W~' cross section is 
very sensitive to the W mass. The sensitivity is largest 
in the region around ^/s = 161 GeV [ at which point 
the statistical uncertainty is given by 

-1/2 



5M^"* « 90 MeV/c^ 



100 pb~ 

Here, e is the efficiency for detecting W bosons. 
e = 0.67 and an integrated luminosity of 100 fb^^, 
finds from Eq. (^) 

SM^"' w 3.5 MeV/c^. 



(8) 

For 
one 

(9) 



Assuming that the efficiency and the integrated 
luminosity can be determined with a precision of Ae = 
0.25% and AC = 0.1%, Mw can be measured with an 
uncertainty of [ H 



provided that the theoretical uncertainty on the 
W^W~ cross section is smaller than about 0.1% in 
the region of interest. 

Presently, the W pair cross section in the threshold 
region is known with an accuracy of about 1.4% [ 
In order to reduce the theoretical uncertainty of 
the cross section to the desired level, the full 0{a) 
electroweak corrections in the threshold region are 
needed. This calculation is extremely difficult. In 
particular, currently no practicable solution of the 
gauge invariance problem associated with finite W 
width effects in loop calculations exists. The existing 
calculations which take into account 0(oi) electroweak 
corrections all ignore non- resonant diagrams [ Q . 

If one (pessimistically) assumes that the theoretical 
uncertainty of the cross section will not improve, the 
uncertainty of the W mass obtained from a threshold 
scan is completely dominated by the theoretical error, 
and the precision of the W mass is limited to [ E^] 



(5M, 



w 



17 MeV/c^ 
24 MeV/c^. 



Act 



X 100% 



(11) 



Using direct reconstruction of W bosons and as- 
suming an integrated luminosity of 500 fb~^ at ^/s = 
500 GeV, one expects a statistical error of SM^""* « 
3.5 MeV/c^ [ Systematic errors are dominated 

by jet resolution effects. Using Z7, Z ^ 2 jet 
events where the photon is lost in the beam pipe for 
calibration, a systematic error ^M^** < 10 MeV/c^ is 
expected to be achieved. The resulting overall preci- 
sion of the W boson mass from direct W reconstruction 
at a Linear Collider operating at an energy well above 
the W pair threshold is 



SMw « 10 MeV/c^ 

9. Theoretical Issues at high y/s 



(12) 



SMw ~ 6 MeV/c^ 



(10) 



Future hadron and lepton collider experiments are 
expected to measure the W boson mass with a preci- 
sion of 6Mw w 10 - 20 MeV/c2. For values of SMw 
smaller than about 40 MeV/c^, the precise definition of 
the W mass and width become important when these 
quantities are extracted. 

In a field theoretical description, finite width effects 
are taken into account in a calculation by resumming 
the imaginary part of the W vacuum polarization. 
This leads to an energy dependent width. However, 
the simple resumming procedure carries the risk of 
breaking gauge invariance. Gauge invariance works 
order by order in perturbation theory. By resumming 
the self energy corrections one only takes into account 
part of the higher order corrections. Apart from being 
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theoretically unacceptable, breaking gauge invariance 
may result in large numerical errors in cross section 
calculations. 

In order to restore gauge invariance, one can adopt 
the strategy of finding the minimal set of Feynman di- 
agrams that is necessary for compensating those terms 
caused by an energy dependent width which violate 
gauge invariance [ This is relatively straightfor- 
ward for a simple process such as qq' W ^ £iy [ Is^ , 



but more tricky for 



w+w- 



4 fermions, in 



particular when higher order corrections are included. 
The so-called complex mass scheme [ |5^, which uses 
a constant, ie. an energy independent width, offers 
a convenient alternative. At LEP II energies, ^/s ~ 
200 GeV, the differences in the e~^e~ — *■ 4 fermions 
cross section using an energy dependent and a constant 
width are small. However, at Linear Collider energies, 
•y/i = 0.5 — 2 TeV, the terms associated with an energy 
dependent width which break gauge invariance lead 
to an overestimation of the cross section by up to a 
factor 3 [HI. 

For qq' ^ W ^ (.v, the parameterizations of the 
W resonance in terms of an energy dependent and a 
constant W width are equivalent. The W resonance 
parameters in the constant width scenario, Mw and 
Tw , and the corresponding quantities, Mw and Fvf, of 
the parameterization using an energy dependent width 
are related by a simple transformation [ p5[ 

-1/2 



Mw 
Tw 



where 7 
constant 



Mw (1 + 7') 

Tw (1+7')" 
= VwlMw- 



1/2 



(13) 
(14) 



The W mass obtained in the 
width scenario thus is about 27 MeV/c^ 
smaller than that extracted using an energy dependent 
width. 

In the past, an energy dependent W width has been 
used in measurements of the W mass at the Tevatron [ 
^ [57| . The Monte Carlo programs available for the 
W mass analysis at LEP II (see Ref. [ ^ for an 
overview) in contrast use a constant W width. Since 
the difference between the W mass obtained using 
a constant and an energy dependent width is of the 
same size or larger than the expected experimental 
uncertainty, it will be important to correct for this 
difference in future measurements. 

10. Conclusions 

The measurements of the W mass and width in 
Run I already represent great experimental achieve- 
ments and contribute significantly to their world av- 
erage determinations. Close inspection of the various 
systematic error sources leads us to believe that a W 
mass measurement in Run II at the 30 MeV/c^ level 



per experiment is achievable, and this compares well 
to the expected uncertainty on the W mass measured 
at LEP II. Each experiment is expected to measure the 
W width to a similar precision with 2 fb~^ of data. 

Alternative methods for determining Mw at the 
Tevatron have been discussed and may turn out to be 
more appropriate in the Run II operating environment 
than the traditional transverse mass fitting approach. 
Determination of the W mass at the LHC will be 
extremely challenging, using detectors that are not op- 
timized for this measurement. A future linear collider 
should do significantly better. Clearly, the W mass 
and width measurements at the Tevatron in Run II 
will remain the best hadron collider determinations of 
these quantities for many years and will compete with 
the best measurements made elsewhere. 
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Measurement of the Forward- Backward Asymmetry in e^e and fi'^fi events 
with D0 in Run II 
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The forward-backward asymmetry of £^£~ events in Run II can yield a measurement of the effective weak mixing angle 
sin^ and can provide a test of the standard model 'y*/Z interference at £'^£~ invariant masses well above the 200 GeV 
center of mass energy of the LEP II collider. The asymmetry at large partonic center of mass energies can also be used to 
study the properties of possible new neutral gauge bosons. We describe an updated study of the forward-backward asymmetry 
and give estimates of the statistical and systematic uncertainties expected in Run II. The prospects for measuring the weak 
mixing angle at the LHC and a linear collider operating at -^s — AIz are also briefly described. 



1. Introduction 

In this note we present an updated study of the 
prospects for measurement of the forward-backward 
asymmetry in pp ^ / Z ^ £^£~ events. This work 
extends our earlier study described in the TeV2000 re- 
port [|l[ in several respects: (i) we include the effects of 
QED corrections; (ii) we include the effects of expected 
Run II D0 detector resolutions and efffciencies; (iii) 
we consider systematic errors in more detail; and (iv) 
we include a simulation of the muon channel process 
pp J* / Z ^ fi+ fi- . 

The forward-backward asymmetry (Apb) in pp —f 
7*/Z events arises from the parton level 

process qq ^ ^* /Z ^ t^l^ . This asymmetry depends 
on the vector and axial-vector couplings of the quarks 
and leptons to the Z boson and is therefore sensitive 
to the effective weak mixing angle sin^ 9w ■ The 
current world average value of sin^ 9w from LEP and 
SLD asymmetry measurements is sin^ 9w = 0.23147 ± 
0.00017 [ H]. As will be seen from our results it will 
be necessary to achieve high luminosity (> 10 fb~^) 
and combine the results from the electron and muon 
channels and the results from D0 and CDF to achieve 
a precision comparable to this. 

The SM tree level prediction [ ^ for ApB as a 
function of s for qq 'y*/Z —>■ e+e" is shown in 
Fig. |l| for u and d quarks. These are the same asymme- 
tries as encountered in the inverse e+e~ annihilation 
reactions. The largest asymmetries occur at parton 
center-of-mass energies of around 70 GeV and above 
110 GeV. At the Z-pole the asymmetry is dominated 
by the couplings of the Z boson and arises from the 
interference of the vector and axial components of 
its coupling. The asymmetry is proportional to the 
deviation of sin^ 9w from j. At large invariant mass. 



the asymmetry is dominated by j* /Z interference and 
is almost constant 0.6), independent of invariant 
mass. 
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Figure 1. The standard model tree level prediction 
of the forward-backward asymmetry as a function of 
e'^e^ invariant mass for uu e+e~ and dd — > e^e^. 



With sufficient statistics in Run II the forward- 
backward asymmetry can be used to measure sin^ 9w , 
which in turn can provide a constraint on the standard 
model complementary to the measurement of the W 
boson mass. The Tevatron also allows measurement 
of the asymmetry at partonic center-of-mass energies 
above the center of mass energy of LEP II. This 
measurement can be used, not only to confirm the 
standard model 'y*/Z interference which dominates 
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in this region, but also to investigate possible new 
phenomena which may alter ApB, such as new neutral 
gauge bosons [ |^ or large extra dimensions [ || . 

CDF have measured the forward-backward asymme- 
try at the Tevatron using e~^e~ pairs in 110 pb""'^ of 
data at = 1.8 TeV [|. They obtain Aps = 0.070± 
0.016 in the mass region 75 GeV < nig+e- < 105 GeV, 
and Afb = 0.43±0. 10 in the region m^+e- > 105 GeV. 
The much larger Run II statistics will enable Afb to 
be measured with an uncertainty reduced by over an 
order of magnitude. 



2. Simulation 

The simulations presented here use the ZGRAD 
Monte Carlo program [0], which includes 0{a) QED 
radiative corrections to the process pp — s- ^*/Z 
t^t". We simulate this process at = 2.0 TeV 
using the MRST parton distributions as our default 
set. Since the radiative corrections are included in 
ZGRAD, wc denote the process of interest by pp 
7*/Z — > ^+£"(7) in the remainder of this paper. The 
ZGRAD program includes real and virtual corrections 
in the initial and final states. 

In our simulations, the effects of detector resolution 
are modeled by smearing the 4-momenta of the parti- 
cles from ZGRAD according to the estimated resolution 
of the Run II D0 detector. We smear the 4-momenta 
of electrons, positrons and photons according to the 
energy resolution ueu of the calorimeters, which have 
been parametrized using constant, sampling and noise 
terms as 



resolution in \/pt is parametrized as: 



= < 



-EM ' 
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Central Calorimeter 
2 

End Calorimeters 
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where we use the parameters relevant for the Run I 
detector, cem = 0.0115, sem = 0.135, and tiem = 
0.43 for the CC, and cem = 0.0100, sem = 0.157, and 
riEM = 0.29 for the EC. With the addition of the 2 T 
solenoidal magnetic field in Run II, only minor changes 
in these parameters are expected. The transverse 
momentum of muons in the Run II detector will be 
measured in the central tracking system, consisting 
of the Central Fiber Tracker (CFT) and the Silicon 
Microstrip Tracker (SMT). The momentum resolution 
of the tracking system has been studied using the 
fast Monte Carlo mcfast. From these studies the 
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where 

L -- 
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Pt\/ L \ sin i 



tan 



<e <ec 

9r<0< 90° 



(2) 



(3) 



Here a = 0.0017 GeV"\ 7 = 0.018, L is the fraction of 
the projection of the track length in the bending plane 
which is measured in the Tracker, and 9c ~ 23° is the 
polar angle beyond which the number of CFT layers 
crossed by a track starts to decrease. The first term 
in Eq. is due to the detector resolution while the 
second term is due to multiple scattering. 

Figure |^ shows the transverse momentum resolution 
as a function of detector pseudorapidity \r]det \ for tracks 
with a pT of 1, 20 and 100 GeV, while Fig. | shows the 
resolution as a function of pr and \r]det\ in the form 
of a contour plot. For central tracks {rj = 0) with 
Pt — 45 GeV, the resolution is (j{pt)/pt = 8%, to be 
compared with the calorimeter energy resolution for 
45 GeV electrons of cr{E)/E = 2.5%. 
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Figure 2. Parametrized transverse momentum 

resolution for the D0 Run II tracking system (dotted 
lines). The solid lines are the results of a simple Monte 
Carlo simulation taken from the D0 SMT Technical 
Design Report [ || . 
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Figure 3. Transverse momentum resolution for 
the central tracking system plotted as a function of 
transverse momentum and pseudorapidity. 



We assume an overall detection efficiency of 75% 
for e+e" events and 65% for n'^fj," events. These 
efficiencies are rough estimates of the effects of trig- 
ger and particle identification efficiencies expected in 
Run II. The results can be updated once more realistic 
numbers for these efficiencies become available. 

The ZGRAD program generates weighted events. Due 
to the occurrence of negative weights, we did not 
attempt to unweight the events. Thus, we work with 
weighted events and properly account for the weights 
in our calculations of the forward-backward asymmetry 
errors. The forward-backward asymmetry is defined by 

A,s = (4) 

where ap and ag are the forward and backward cross 
sections, defined by 



da 



d{cose* 
da 



-d{cosd*) 



and 9* is the angle of the lepton in the CoUins-Soper 
frame [ ^ . 

The statistical error on Aps is given by 



SApB = 2 



[crp + aBY ^ ' 

where Sap, 5a b are the uncertainties in the forward 
and backward cross sections. For unweighted events. 



this simplifies to 
2 



5A 



FB 



Np + Nb \ Np + N, 



NpNb 



(7) 



where Np, Nb are the numbers of forward and 
backward events. However, ZGRAD generates weighted 
events and, therefore, we use Eq. (||) where Sap 
and SaB, are calculated using the appropriate event 
weights. 

The selection cuts used in our study are summarized 
in Table In the electron channel we require one of the 
electrons to be in the CC (|??det| < 1-0), while the other 
electron may be in the CC or in the EC (|77ciet| < 1-0 
or 1.5 < \T]det\ < 2.5). 

In the muon channel we require both muons to be 
within \r]det\ < f.7. In Run II the muon coverage is 
expected to extend up to \r]det\ — 2.0. We chose to limit 
the muon acceptance to \r]det\ = 1-7 since Monte Carlo 
events were already generated with this restriction 
and large CPU time would have been required to 
re-generate the events. 

We account in our simulation for the granularity 
of the D0 calorimeter. If the photon is very close 
to the electron its energy will be merged with that 
of the electron cluster. Thus, in the simulation we 
combine the photon and electron 4-momenta to form 
an effective electron 4-momentum if the photon is 

within AR. 



1.67 — \/ ^Ve-y 



A02 



falls within 0.2 < AR, 



'67 



< 0.2. If the photon 

< 0.4, we reject the event if 
E^/ {Ee + Ej) > 0.15, since the event will not pass the 
standard isolation criterion imposed on electrons. 

If a photon is very close to a muon and it deposits 
sufficient energy in the calorimeter close to the muon 
track, the energy deposition in the calorimeter will not 
be consistent with the passage of a minimum ionizing 
muon. Therefore, in the simulation we reject events if 
AR^,-y < 0.2 and > 2 GeV. 



Table 1 

Selection criteria for e+e~ and fJ-^ fi~ events. 







e+e 






Selection cut 


ei 


62 


Ml 




PT (GeV) 


> 25 


> 25 


> 20 


> 15 


\r]det\ 


< 1.0 


< 1.0 


< 1.7 


< 1.7 






or 1.5 - 2.5 






nii+i- (GeV) 


> 40 


> 40 


> 40 


> 40 



25 



3. Results 

The t' invariant mass distributions for pp 
-f*/Z i+i~{-f) at Vs = 2.0 TeV from the 

ZGRAD simulations, using the MRST parton distri- 
bution functions are shown in Fig. ^. The thin hne 
shows da/dmi+g- without any kinematic cuts apphed 
and with no detector acceptance or resolution effects 
included. fn order to obtain sufficient statistical 
precision a large number of events were generated in 
multiple runs covering overlapping regions of m^+i-. 
The thick line shows da/dm^+g- after kinematic cuts 
and detector effects are included. The error bars 
represent the statistical errors only, calculated from 
Eq. (H), assuming an integrated luminosity of f fb^^. 

Fig. ^ shows the forward-backward asymmetry as a 
function of m£+£-. The solid line shows Aps without 
any kinematic cuts applied and with no detector accep- 
tance or resolution effects included. The solid points 
show ApB after kinematic cuts and detector effects are 
included. The error bars represent the statistical errors 
only, calculated from Eq. (|6|), assuming an integrated 
luminosity of 10 fb^^. As can be seen, detector 
resolution and acceptance effects significantly alter the 
shape of the Aps vs. mg+i- curve, especially at low 
di-lepton invariant masses. In this region, the effect 
of CC/EC acceptance increases ApB, while restricting 
events to be in the central region decreases 
the asymmetry. This is also true of e^e~ events if 
only CC/CC events are considered. In the vicinity 
of the Z-pole the energy resolution is better than the 
Pt resolution, and hence the A'^p^ is altered less 

than A^g . In these plots the ApB shown is the 
reconstructed ApB without corrections for acceptance 
or resolution effects. 

In order to obtain a measurement of the weak mixing 
angle we assume the relationship 
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FB 



b sin 



(8) 



so that the statistical error on sin 0w is given by 
5ApB 



S sin 



yw — 



(9) 



The quantity b is determined by varying sin 6w in the 
Monte Carlo simulations. Since ApB is determined 
over a finite range of di-lcpton invariant mass, we 
have investigated the effect of the lower and upper 
m£+£- cuts on 6sin^9y[^. The optimal precision is 
obtained for 75 GeV < nif+g- < 105 GeV, i.e. a 
mass window encompassing the Z-pole. This is to be 
expected because the sensitivity b is maximal at the 
Z-pole and this region is where the cross section peaks 
and hence the statistical error is smallest. Thus, the 
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Figure 4. Invariant mass mi+i- distributions for (a) 
e^e~ events and (b) fJ-'^ fJ.~ events. The thin line is the 
distribution obtained with no cuts or detector effects 
applied and the thick line is the resulting distribution 
after selection cuts and detector effects are included. 



ApB values and errors presented in the remainder of 
this paper are all obtained with a di-lepton invariant 
mass cut of 75 GeV < m£+^- < 105 GeV. Table || 
shows the resulting statistical uncertainties obtained 
from the electron and muon channels. The e^e~ and 
fj-^ IJ-^ channels yield similar uncertainties on Aps and 
sin^ 9w ■ In both channels the effect of the selection 
cuts is to reduce the sensitivity b from about 3.5 to 
about 2.8. 

The effects of NLO QCD corrections to the process 
PP ~* £^£^ are not included in ZGRAD, so 

we estimate these using the 0{as) event generator 
described in [ Q . Using the same method as described 
in Section ^, we calculate the change in the forward- 
backward asymmetry and the shift in sensitivity due 
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Table 2 

Uncertainties on Aps and s'm^ 9w in the invariant mass range 75 GeV < m^+i- < 105 GeV for an integrated 
luminosity of 10 fb~^. Also shown are the assumed event detection efficiency, the number of events passing all the 
cuts A^obs, the forward-backward asymmetry ApB, and the sensitivity h. 



Process 
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Efficiency 


Nobs 


Afb 


s A Stat 
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S sin^ 9w 


Z e+e- 


no 


100% 


1.78 X 10" 


0.0551 


0.0008 


3.43 


0.0002 


Z e+e- 


yes 


75% 


3.82 X 10^ 


0.0515 


0.0014 


2.78 


0.0005 
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no 


100% 


1.87 X 10^ 


0.0534 


0.0007 


3.51 


0.0002 


Z — !■ fi+^~ 


yes 


65% 


5.67 X 10'"^ 


0.0420 


0.0011 


2.62 


0.0004 



to NLO QCD corrections. Thus, we write 
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Figure 5. Forward-backward asymmetry Afb vs. 
di-lepton invariant mass for (a) e+e~ events and 
(b) events. The solid line is the distribution 

obtained with no cuts or detector effects applied and 
the points arc the resulting distribution after selection 
cuts and detector effects are included. The error bars 
represent the statistical errors for a data sample of 
10 fb-i. 



A6 



= A 



FB 



-A 



LO 
FB 
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where LO denotes the leading-order quantities. For 
events generated including detector effects we find the 
shift in Afb to be negligible for e+e- events and 
— 13% for iJi+ jJT events. The shift in sensitivity is 
A6/6^o « -3.4% for e+e- events and « -25% for 
(/i+/i-) events. Thus, NLO QCD effects decrease the 
sensitivity to sin^ 9w by 3.4% in the e+e- channel and 
25% in the /i+/J.~ channel. 

4. Systematic Uncertainties 

4.1. Parton Distribution Functions 

Since the vector and axial couplings of the u and 
d quarks to the Z boson are different, the lepton 
forward-backward asymmetry is expected to depend 
on the ratio of the u to d quark parton distribution 
functions. Thus, the choice of the parton distribution 
functions (PDF's) will affect the measured lepton 
forward-backward asymmetry. 

We have run simulations with six PDF's from the 
MRS [ |l§ and CTEQ [ |ll| sets to study the effect of 
the PDF's on the asymmetry. Fig. |6| shows the e+e~ 
and /i+/i- asymmetries and their statistical errors for 
each PDF. 

The largest deviation from the MRST value for Afb 
is 0.0018 for the e+e~ channel and 0.0015 for the fJ,+^J,- 
channel. While these numbers are of the same order 
as the statistical error expected on Afb for 10 fb-^, 
we expect that in Run II our knowledge of the PDF's 
will improve considerably, e.g. from the constraints 
imposed by the Run II W asymmetry measurements. 
Thus, we expect a significantly decreased systematic 
error due to the uncertainty in the PDF's which 
will likely render it insignificant compared with the 
statistical error in the measurement. For example, if 
the PDF uncertainty scales as 1/ ViV, the uncertainty 
in Afb would be SAfb ~ 0.00018 (0.00015) for an 
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Variation of Forward-Baci<ward Asymmetry with PDF 
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Figure 6. Variation of the forward-backward asym- 
metry for e~^e~ and fJ.^ IJ-^ events, including selection 
cuts and detector effects. 
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4.2. Energy Scale Calibration Uncertainties 

The energy scale uncertainty, or the uncertainty in 
mapping the calorimeter response to the true electron 
energy, affects the forward-backward asymmetry by 
causing a shift in the £~^£~ invariant mass range over 
which we integrate Aps- The effect is significant in 
the Z-pole region, but at large invariant masses ApB 
is essentially constant and the energy scale uncertainty 
does not play a role. For electrons, the measured 
energy can be related to the true energy by 



E„ 



= aEtr 



(12) 



where the scale factor a and offset S are determined 
by calibration of the calorimeters. In Run I D0 deter- 
mined a = 0.9533 ± 0.0008 and S = -0.16tg;^i^ GeV. 

Assuming these uncertainties in a and S, we find a 
systematic error of SApB = 0.0002 due to the overall 
energy scale uncertainty. 

4.3. Uncertainty due to Backgrounds 

Backgrounds are not included in the simulations 
above, but we can estimate the uncertainty due to 
backgrounds as follows. If we assume that the fraction 
of observed events which are due to backgrounds is 
a ± (5a, then the uncertainty in the forward-backward 
asymmetry will be 



SApB = 2V2 



NpN, 



B 



{Np + Ne 



(13) 



where Np, Nb are the numbers of background- 
subtracted forward and backward events, and we have 
assumed that the background events are symmetric in 
cos9* . If we assume that the error in the background 
fraction 5a scales as the inverse of the integrated 
luminosity, we can extrapolate from the uncertainties 
in the Run I CDF and D0 data samples to estimate the 
error. The Run I uncertainties were 5a « 0.1 — 1.0%. 
Thus, for an integrated luminosity of 10 fb~^, we 
obtain an uncertainty on ApB of 5ApB = 0.00014. 

4.4. Summary of Uncertainties in sin^ 6w 
Table ^ summarizes the statistical and individ- 
ual systematic uncertainty estimates expected with 
10 fb~^ of data. We estimated the uncertainties due 
to electron energy resolution and muon transverse 
momentum resolution to be negligible. 

4.5. Conclusions 

The measurement of the forward-backward asymme- 
try in e+e~ and /i^/i^ events in Run II provides a 
means to test the standard model "f* /Z interference at 
£~^i~ invariant masses well above the center of mass 
energy of the LEP II collider. The estimated D0 
precision on ApB achievable with 10 fb~^ integrated 
luminosity is shown in Fig. |^. 

In the vicinity of the Z-po\e this measurement can 
also be used to determine the effective weak mixing 
angle sin^ 9w . The optimal precision on sin^ 9w is 
obtained for 75 GeV < mg+i^ < 105 GeV, i.e. a 
mass window encompassing the Z-po\e. This is to 
be expected because the sensitivity b is maximal at 
the Z-pole and this region is where the cross section 
peaks and hence the statistical error is smallest. For 
10 fb^^ we estimate that the total error on sin^ 9w 
will be 0.0005 in the electron channel and 0.0006 in 
the muon channel, assuming that systematic errors 
scale as the inverse of the square root of the integrated 
luminosity. One would expect similar precision from 
CDF, and combining the results of the two experiments 
in both channels the overall uncertainty would be 
5sm^ 9w ~ 0.00028. Therefore, if integrated luminosi- 
ties in excess of 10 fb~^ can be achieved in Run II, 
it appears that the determination of sin^ 9w will have 
comparable precision to the current world average of 
the measurements from LEP and SLD. 



5. Measuring ApB at the LHC 

At the LHC, the Z £^£~ cross section is approxi- 
mately a factor 7 larger than at the Tevatron. However, 
the measurement of the forward backward asymmetry 
is complicated by several factors. In pp collisions, the 
quark direction in the initial state has to be extracted 
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Table 3 

Summary of uncertainties on Aps and s'm^ 9^ in the invariant mass range 75 GeV < m£+£- < 105 GeV for an 
integrated luminosity of 10 fb^^. The effects of 0{a) QED corrections and NLO QCD corrections have been taken 
into account. 





e+ 


e 






Error source 


SAfb 


S sin 6w 


SAfb 


S sin^ 9w 


Statistical 


0.0014 


0.0005 


0.0011 


0.0006 


Systematics 


0.0003 


0.0001 


0.0002 


0.00010 


PDF 


0.00018 


0.00007 


0.00015 


0.00008 


EM scale 


0.0002 


0.00007 






Backgrounds 


0.00014 


0.00005 


0.00014 


0.00007 


Total uncertainty 


0.0014 


0.0005 


0.0011 


0.0006 



from the boost direction of the £^£^ system with 
respect to the beam axis. At LHC energies, the sea-sea 
quark flux is much larger than at the Tevatron. As a 
result, the probability, fq, that the quark direction and 
the boost direction of the di-lepton system coincide 
is rather small. The forward backward asymmetry 
is therefore smaller than at the Tevatron, and the 
sensitivity to sin^ 9w at the LHC with 100 fb~^ per 
lepton channel and experiment [ 0, |l^ is similar to that 
estimated for the Tevatron with 10 fb^^ (see Sec. 1.5). 

Restricting the Aps measurement to events which 
satisfy \y{i+£-)\ > 1 in addition to the \y{i)\ < 2.5 
cut improves the significance of the measurement by 
about a factor 1.5. Events with a large di-lepton 
rapidity originate from collisions where at least one 
of the partons carries a large fraction x of the proton 
momentum. Since valence quarks dominate at high 
values of x, a cut on y{£~^£~) increases fq and thus the 
asymmetry. However, the gain due to the larger asym- 
metry is partially cancelled by the loss of statistics, 
leaving a modest improvement only. 

In order to achieve a precision better than the 
current value of S sin^ 9w = 1.7 X 10-4 [|l|l, it wiU be 
necessary to detect one of the leptons in the rapidity 
range up to \y{()\ < 5 at the LHC. If this can be 
done, one expects that the weak mixing angle can be 
determined with a precision of 



d sin^ 9 



w 



1.4 X 10" 



(14) 



per lepton channel and experiment for an integrated 
luminosity of 100 fb^^. In order to reach the precision 
given in Eq. (14), a jet rejection factor of 10 — 100 
has to be achieved in the forward rapidity region 
2.5 < \y{£)\ < 5, and the lepton acceptance times the 
reconstruction efficiency as a function of y{i) has to be 
known to 0.1% or better [ H). 

For comparison, at a Linear Collider operating at 



y/s = Mz with a luminosity of a few x 10^'^ cm"^ s~^, 
it is expected that the weak mixing angle be deter- 
mined with a precision of about Ssin^ 9w = 1 x 10"^ [ 

HI. 
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Global Fits to Electroweak Data Using GAPP 
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At Run II of the Tevatron it will be possible to measure the W boson mass with a relative precision of about 2 x 10~*, 
which will eventually represent the best measured observable beyond the input parameters of the SM. Proper interpretation 
of such an ultra-high precision measurement, either within the SM or beyond, requires the meticulous implementation and 
control of higher order radiative corrections. The FORTRAN package GAPP, described here, is specifically designed to meet 
this need and to ensure the highest possible degrees of accuracy, reliability, adaptability, and efficiency. 



1. PRECISION TESTS 

Precision analysis of electroweak interactions follows 
three major objectives: high precision tests of the 
SM; the determination of its fundamental parameters; 
and studies of indications and constraints of possible 
new physics beyond the SM, such as supersymmetry 
or new gauge bosons. Currently, the experimental 
information comes from the very high precision Z 
boson measurements at LEP and the SLC, direct mass 
measurements and constraints from the Tevatron and 
LEP II, and low energy precision experiments, such 
as in atomic parity violation, p scattering, and rare 
decays. These measurements are compared with the 
predictions of the SM and its extensions. The level of 
precision is generally very high. Besides the need for 
high-order loop calculations, it is important to utilize 
efficient renormalization schemes and scales to ensure 
sufficient convergence of the perturbative expansions. 

The tasks involved called for the creation of a 
special purpose FORTRAN package, GAPP, short 
for the Global Analysis of Particle Properties [ |^. 
It is mainly devoted to the calculation of pseudo- 
observables, i.e., observables appropriately idealized 
from the experimental reality. The reduction of 
raw data to pseudo-observables is performed by the 
experimenters with available packages (e.g., ZFITTER 
for Z pole physics). For cross section and asymmetry 
measurements at LEP II (not implemented in the 
current version, GAPP_99.7), however, this reduction 
is not optimal and convoluted expressions should be 
used instead. GAPP attempts to gather all available 
theoretical and experimental information; it allows the 
addition of extra parameters describing new physics; it 
treats all relevant SM inputs as global fit parameters; 
and it can easily be updated with new calculations, 
data, observables, or fit parameters. For clarity and 
speed it avoids numerical integrations throughout. It 
is based on the modified minimal subtraction (MS) 



scheme which demonstrably avoids large expansion 
coefficients. 

GAPP is endowed with the option to constrain 
nonstandard contributions to the oblique parameters 
defined to affect only the gauge boson self-energies [ |j 
(e.g. S, T, and U); specific anomalous Z couplings; the 
number of active neutrinos (with standard couplings to 
the Z boson); and the masses, mixings, and coupling 
strengths of extra Z bosons appearing in models 
of new physics. With view on the importance of 
supersymmetric extensions of the SM on one hand, and 
upcoming experiments on the other, I also included 
the 6 — !■ S7 transition amplitude, and intend to 
add the muon anomalous magnetic moment. In the 
latter case, there are theoretical uncertainties from 
hadronic contributions which arc partially correlated 
with the renormalization group (RG) evolutions of the 
QED coupling and the weak mixing angle. These 
correlations will be partially taken into account by 
including heavy quark effects in analytical form; see 
Ref. [ ^ for a first step in this direction. By comparing 
this scheme with more conventional ones, it will also be 
possible to isolate a QCD sum rule and to rigorously 
determine the charm and bottom quark MS masses. 
The and rhb, with high precision. 

2. GAPP 

2.1. Basic structure 

In the default running mode of the current version, 
GAPP_99.7, a fit is performed to 41 observables, out 
of which 26 are from Z pole measurements at LEP 
and the SLC. The Fermi constant, Gp (from the muon 
lifetime), the electromagnetic fine structure constant, 
a (from the quantum Hall effect), and the light fermion 
masses are treated as fixed inputs. The exception is 
rhc which strongly affects the RG running of a{^) 
for fi > rhc- I therefore treat rhc as a fit parameter 

^Quantities defined in the MS scheme are denoted by a caret. 
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and include an external constraint with an enhanced 
error to absorb hadronic threshold uncertainties of 
other quark flavors, as well as theoretical uncertainties 
from the application of perturbative QCD at relatively 
low energies. Other fit parameters are the Z boson 
mass, Mz, the Higgs boson mass, Mjy, the top quark 
mass, mt, and the strong coupling constant, a^, so 
that there are 37 effective degrees of freedom. Given 
current precisions, Mz may alternatively be treated as 
an additional fixed input. 

The file fit.f basically contains a simple call to 
the minimization program MINUIT [ H (from the 
CERN program library) which is currently used in data 
driven mode (see smfit.dat). It in turn calls the 
core subroutine fen and the x^-function chi2, both 
contained in chi2 . f . Subroutine fen defines constants 
and flags; initializes parts of the one-loop package FF [ 
^ ^; and makes the final call to subroutine values 
in main . f which drives the output (written to file 
smfit.out). In chi2 the user actively changes and 
updates the data for the central values, errors, and 
correlation coefficients of the observables, and includes 
or excludes individual contributions to (right after 
the initialization, clii2 = O.dO). To each observable 
(as defined at the beginning of chi2) corresponds an 
entry in each of the fields value, error, smval, and 
pull, containing the central observed value, the total 
(experimental and theoretical) error, the calculated fit 
value, and the standard deviation, respectively. The 
function chi2 also contains calls to various other sub- 
routines where the actual observable calculations take 
place. These are detailed in the following subsections. 

Another entry to GAPP is provided through mh.f 
which computes the probability distribution function 
of Mh. The probability distribution function is the 
quantity of interest within Bayesian data analysis (as 
opposed to point estimates frequently used in the con- 
text of classical methods) , and defined as the product 
of a prior density and the hkelihood, L ~ exp(— x^/2). 
If one chooses to disregard any further information on 
Mh (such as from triviality considerations or direct 
searches) one needs a non-informative prior. It is 
recommended to choose a flat prior in a variable 
defined on the whole real axis, which in the case of 
Mh is achieved by an equidistant scan over logMjy. 
An informative prior is obtained by activating one of 
the approximate Higgs exclusion curves from LEP II 
near the end of chi2.f. These curves affect values 
of Mh even larger than the corresponding quoted 
95% GL lower limit and includes an extrapolation to 
the kinematic limit; notice that this corresponds to a 
conservative treatment of the upper Mh limit. 

Contour plots can be obtained using the routine 
mncontours from MINUIT. For the cases this fails, 



some simpler and slower but more robust contour 
programs are also included in GAPP, but these have 
to be adapted by the user to the case at hand. 

2.2. a, sin^ 9w, Mw 

At the core of present day electroweak analyses is 
the interdependence between Gp, Mz, the W boson 
mass, Mn/, and the weak mixing angle, sin^^vv- In 
the MS scheme it can be written as [ 0: [s) j 



A' 



Afvv)' 



-2 -2 
S C 



Af|(l-Afz)' 



where. 



A = 



-,1/2 



37.2805(2) GeV, 



(1) 



(2) 



is the MS mixing angle, 



Afiv 



and, 



-A. 



n 



WW 



(M^)-nw(o) 



V + B, (3) 



Afz = Afvv -f (1 - Afvy) 



nzz(M|) 



Ml 



(4) 



collect 



the 



radiative corrections 
computed in sin2th. f . The II indicate MS subtracted 
self-energies, and V + B denote the vertex and box 
contributions to decay. Although these relations 
involve the MS gauge couplings they employ on-shell 
gauge boson masses, absorbing a large class of radiative 
corrections [ ^ . 

lS.fw and Af^ are both dominated by the contribu- 
tion l^~f{Mz) which is familiar from the RG running 
of the electromagnetic coupling, 

Q. 

aiu) = , (5) 

l-fA,(/x)' 

and computed in alfahat.f up to four-loop 0{aa^^). 
Gontributions from c and b quarks are calculated using 
an unsubtracted dispersion relation [ ^ . If is equal to 
the mass of a quark, three-loop matching is performed 
and the definition of a changes accordingly. Pure QED 
effects are included up to next-to- leading order (NLO) 
while higher orders are negligible. Precise results can 
be obtained for /i < 2m.,: and /i > m,c. 

Besides full one-loop electroweak corrections, /S.fy/ 
and Af^ include enhanced two-loop contributions 
of 0{a^n4) [0 (implemented using the ana- 
lytic expressions of Ref. [ |l|) and 0{a'^mf) [ |l2| 
(available as expansions in small and large Mh); 
mixed electroweak/ QGD corrections of 0{aas) [ |l3| 
and ©(aa^m?) [ |l|; the analogous mixed elec- 
troweak/QED corrections of 0{a^); and fermion mass 
corrections also including the leading gluonic and 
photonic corrections. 
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2.3. Z decay widths and asymmetries 

The partial width for Z f f decays is given by, 



_ N[;Mza 



Pf\ [l-4|Q/|Re(%)s2 + 8g 



1 + (Jqed 



cNS 

"qcd 



^QCD 



(6) 



/ 
c 
and 



N. 



is the color factor, Qf is the fermion charge, 
Pf and kf are form factors which differ from 



unity through one-loop electroweak corrections [ 
and are computed in rho . f and kappa . f , respectively. 
For f ^ b there are no corrections of 0{a^mj) and 
contributions of 0{a'^m^) to kf [ |l^ and pf [ || are 
very small and presently neglected. On the other hand, 
vertex corrections of 0{aas) [ |l^ are important and 
shift the extracted as by ~ 0.0007. 

The Z ^ bb vertex receives extra corrections due 
to heavy top quark loops. They are large and have 
been implemented in bvertex.f based on Ref. [ [l^. 
0{a'^m*) corrections [ are included, as well, 
while those of 0{a'^mf) are presently unknown. The 
leading QCD effects oiO{aasm^) and all sublead- 
ing 0{aas) corrections [ ^ are incorporated into pb 
and Kb, but not the O^aa^rn^) contribution which is 
presently available only for nonsinglet diagrams [ ^ . 

In Eq. d), (5qed are the 0{a) and 0{a^) QED 
corrections. i5qcd ^^'^ universal QCD corrections 
up to 0{al) which include quark mass dependent 
contributions due to double-bubble type diagrams [ 
p2| |2^ . <5qcd are the singlet contributions to the 
axial- vector and vector partial widths which start, 
respectively, at 0{al) and 0{a^), and induce relatively 
large family universal but flavor non-universal rrit 
effects [ |24|. The corrections appearing in the 
second line of Eq. are evaluated in leplOO.f . 

The dominant massless contribution to Sqq^ can 
be obtained by analytical continuation of the Adler 
D-function, which (in the MS scheme) has a very well 
behaved perturbative expansion ^ 1 + X]i=o ^i'^s^^ 
Qs — a.s{Mz)/Tr (see the Appendix for details). The 
process of analytical continuation from the Euclidean 
to the physical region induces further terms which 
are proportional to /3-function coefficients, enhanced 
by powers of tt^, and start at 0{a^g). Fortunately, 
these terms [ ^ involve only known coefficients up to 
©(ofg), and the only unknown coefficient in ©(a^) is 
proportional to the four-loop Adler function coefficient, 
da. In the massless approximation. 



"QCD 



as + 1.4092a2 - (0.681 + 12.086) 
(da - 89.19) ai + {di + 79.7) + 
(ds - 121(^3-^3316) al, 



(7) 



and terms of order a], ^ lO""'^" are clearly negligible. 
Notice, that the C(as) term effectively reduces the 
sensitivity to d^ by about 18%. Eq. (|^) amounts to a 
reorganization of the perturbative series in terms of the 
di times some function of as ; a similar idea is routinely 
applied to the perturbative QCD contribution to t 
decays [ ^ . 

Final state fermion mass effects [ 22, |2^ of 0{rnj) 
(and 0{mf) for b quarks) are best evaluated by 
expanding in Th'^{Mz) thus avoiding large logarithms 
in the quark masses. The singlet contribution of 
0{alinl) is also included. 

The dominant theoretical uncertainty in the Z line- 
shape determination of as originates from the massless 
quark contribution, and amounts to about ±0.0004 
as estimated in the Appendix. There are several 
further uncertainties, all of O{10~^): from the 0(as) 
heavy top quark contribution to the axial-vector part 
of ^QCD' from the missing 0{aalm^) and 0{a^m^) 
contributions to the Zfcfe-vertex; from further non- 
enhanced but cohering 0(aQ;g)-vertex corrections; and 
from possible contributions of non-perturbative origin. 
The total theory uncertainty is therefore. 



Aas(Afz) = ±0.0005, 



(8) 



which can be neglected compared to the current ex- 
perimental error. If rhi, is kept fixed in a fit, then its 
parametric error would add an uncertainty of ±0.0002, 
but this would not change the total uncertainty (||). 

Polarization asymmetries are (in some cases up to 
a trivial factor 3/4 or a sign) given by the asymmetry 
parameters, 



A 



f 



l~4\Qf\Reikf)s^ 



1 - i\Qf\Re{kf)P + 8Q}S^kf\^ 



and the forward-backward asymmetries by. 



AfbU) = -^AeAf. 



(9) 



(10) 



The hadronic charge asymmetry, Qfb, is the linear 
combination. 



q—d^s,b 



RqApBiq)- RqA-Fsiq), 



(11) 



q—u^c 



and the hadronic peak cross section, (Thad, is stored in 
sigmah, and defined by. 



127rr, 



fhad 



Fhad 



Af|r| 



(12) 



Widths and asymmetries are stored in the fields 
gamma(f), alr(f), and afb(f). The fermion index, 
f, and the partial width ratios, R(f), are defined in 
Table 0. 
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Table 1 

Some of the variables used in leplOO.f. Finv and 
Thad are the invisible and hadronic decays widths, 
respectively. 
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gamma (0) = Tiny 


alr(O) = 


1 


1 


e 


R(l) = Thad/r^+e- 


alr(l) = 




2 




R(2) = rhad/r^+^- 


alr(2) = 
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R(3) = Fhad/rT + T-- 


alr(3) = 


Ar 
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R(4) = r^/Fhad 
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R(5) = Fcc/Fhad 


alr(5) = 


Ac 


6 
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7 


d 


R(7) = r<irf-/rhad 
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R(8) = Fss/Fhad 


alr(8) = 


As 
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h 


R(9) = Tfcjrhad 


alr(9) = 


Ab 


10 


had 


gamma (10) Thad 


afb(lO) = 


- Qfb 


11 


all 


gcunmaCll) = Tz 







2.4. Fermion masses 

I use MS masses as far as QCD is concerned, 
but retain on-shell masses for QED since renormalon 
effects are unimportant in this case. This results in 
a hybrid definition for quarks. Accordingly, the RG 
running of the masses to scales /i ^ rhq uses pure 
QCD anomalous dimensions. The running masses 
correspond to the functions msrun (/i) , mcrun (^) , etc. 
which are calculated in masses. f to three-loop order. 
Anomalous dimensions are also available at four-loop 
order [ but can safely be neglected. Also needed 
is the RG evolution of as which is implemented to 
four-loop precision in alfas.f. 

I avoid pole masses for the five light quarks through- 
out. Due to renormalon effects, these can be deter- 
mined only up to ©(Aqcd) and would therefore induce 
an irreducible uncertainty of about 0.5 GeV. In fact, 
perturbative expansions involving the pole mass show 
unsatisfactory convergence. In contrast, the MS mass 
is a short distance mass which can, in principle, be 
determined to arbitrary precision, and perturbative 
expansions are well behaved with coefficients of order 
unity (times group theoretical factors which grow only 
geometrically). Note, however, that the coefficients of 
expansions involving large powers of the mass, m", are 
rather expected to be of 0{n). This applies, e.g., to 
decays of heavy quarks {n — 5) and to higher orders in 
light quark mass expansions. 

The top quark pole mass enters the analysis 
when the results on mt from on-shell produced top 
quarks at the Tevatron are included. In subroutine 
polemasses (nf ,mpole) rhqirhq) is converted to the 
quark pole mass, mpole, using the two- loop pertur- 
bative relation from Ref . [ |3^ . The exact three- loop 
result [ pl| has been approximated (for mt) by em- 



ploying the BLM [ 32 scale for the conversion. Since 
the pole mass is involved it is not surprising that 
the coefficients are growing rapidly. The third order 
contribution is 31%, 75%, and 145% of the second 
order for rrit (nf = 6), uib (nf = 5), and rric (nf = 4), 
respectively. I take the three-loop contribution to 
the top quark pole mass of about 0.5 GeV as the 
theoretical uncertainty, but this is currently negligible 
relative to the experimental error. At a high energy 
lepton collider it will be possible to extract the MS 
top quark mass directly and to abandon quark pole 
masses altogether. 

2.5. ly scattering 

The ratios of neutral-to-charged current cross sec- 
tions, 



Ri, — 



^NC 

,cc ■ 



'vN 
^CC ' 



(13) 



have been measured precisely in deep inelastic v 
(D) hadron scattering (DIS) at CERN (CDHS and 
CHARM) and Fermilab (CCFR). The most precise 
result was obtained by the NuTeV Collaboration at 
Fermilab who determined the Paschos-Wolfenstein ra- 
tio, 



NC _ NC 

JL — r^r^ r~^(~< -1^1/ ' -^y -j 



^cc 

'vN 



(14) 



'vN 



with r — <J^N /'^vN- Results on R^ are frequently 
quoted in terms of the on-shell weak mixing angle 
(or Mw) as this incidentally gives a fair description 
of the dependences on mt and Mh- One can write 
approximately. 



K = gi 

where, 

2 1 



9Rr, 



Ru = gl 



^, R- 

r 



= gl-gl, (15) 



9 



9r = -^sva^Ow (16) 



However, the study of new physics requires the im- 
plementation of the actual linear combinations of 
effective four-Fermi operator coefficients, eL,R{u) and 
^L,R.{d), which have been measured. With the ap- 
propriate value for the average momentum trans- 
fer, q2, as input, these are computed in the sub- 
routines nuh(q2,epsu_L,epsd_L,epsu_R,epsd_R) (ac- 
cording to Ref. [ ^), nuhnutev, nuhccfr, and 
nuhcdhs, all contained in file dis.f. Note, that 
the CHARM results have been adjusted to CDHS 
conditions [ |3^ . While the experimental correlations 
between the various DIS experiments are believed to 
be negligible, large correlations are introduced by the 
physics model through charm mass threshold effects, 
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quark sea effects, radiative corrections, etc. I con- 
structed the matrix of correlation coefficients using the 
analysis in Ref. [ pil , 



It is given by [ 



R- 



Ru Rv Ry Ry Ru Ri/ 



1.00 
0.10 
0.10 
0.10 
0.00 
0.00 
\ 0.00 



0.10 
1.00 
0.40 
0.40 
0.10 
0.10 
0.10 



0.10 
0.40 
1.00 
0.40 
0.10 
0.10 
0.10 



0.10 
0.40 
0.40 
1.00 
0.10 
0.10 
0.10 



0.00 
0.10 
0.10 
0.10 
1.00 
0.15 
0.15 



0.00 
0.10 
0.10 
0.10 
0.15 
1.00 
0.15 



0.00 \ 
0.10 
0.10 
0.10 
0.15 
0.15 
1.00 j 



.(17) 



The effective vector and axial-vector couplings, 
and , from elastic ve scattering are calculated in 
subroutine nue(q2,gvnue,ganue) in file nue.f. The 
momentum transfer, q2, is currently set to zero [ 
Needed is the low energy p parameter, rhonc, 
which describes radiative corrections to the neutral- 
to-charged current interaction strengths. Together 
with sin2t0 (described below) it is computed in file 
lowenergy . f . 

2.6. Low energy observables 

The weak atomic charge, Qw, from atomic parity 
violation and fixed target ep scattering is computed 
in subroutine apv(Qw,Z,AA,Clu,Cld,C2u,C2d) where 
Z and AA are, respectively, the atomic number and 
weight. Also returned are the coefficients from lepton- 
quark effective four-Fermi interactions which are cal- 
culated according to [ |7) . 

These observables are sensitive to the low energy 
mixing angle, sin2t0, which defines the electroweak 
counterpart to the fine structure constant and is similar 
to the one introduced in Ref. [ |^ . There is significant 
correlation between the hadronic uncertainties from 
the RG evolutions of a and the weak mixing angle. 
Presently, this correlation is ignored, but with the 
recent progress in atomic parity violation experiments 
it should be accounted for in the future. 

An additional source of hadronic uncertainty is in- 
troduced by 7.Z-box diagrams which are unsuppressed 
at low energies. At present, this uncertainty can be 
neglected relative to the experimental precision. 

Besides apv, the file pnc . f contains in addition 
the subroutine moller for the anticipated polarized 
fixed target M0ller scattering experiment at SLAG. 
Radiative corrections are included following Ref. [ |3^ . 

2.7. 6 ^ S7 

Subroutine bsgamma returns the decay ratio. 
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(18) 



TT 



Vr 



cb 
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(19) 



where | Vj* Vtf,/Vc6p = 0.950 is a combination of 
Gabbibo-Kobayashi-Maskawa matrix elements and S 
is the Sudakov factor [ 5np and 6ew ^re non- 
perturbative and NLO electroweak [ Q corrections, 
both for the b —i- sj and the semileptonic (6 — ^ cei') 
decay rates. 



47r 



V), 



(20) 



is called the reduced amplitude for the process b S7, 
and is given in terms of the Wilson coefficient C7 
at NLO. Cr and the other d appearing below are 
effective Wilson coefficients with NLO RG evolution [ 
Q from the weak scale to /i = rhb understood. The 
NLO matching conditions at the weak scale have been 
calculated in Ref. [ Q . D includes the virtual gluon 
corrections. 



V^r2C°+rTC^+rsCl 



(21) 



so that it squares to a positive definite branching 
fraction. On the other hand, the amplitude for gluon 
Bremsstrahlung (6 — > sjg) , 



A 



Qa(™ti) 



[CKCl f2sil) + C?/27(l) + C2"/22(l))- 



(22) 



is added linearly to the cross section. The Wilson 
coefficient C2 is defined as in Ref. [ ^ . It enters only 
at NLO, is significantly larger than Cj, and dominates 
the NLO contributions. The parameter < 5 < 1 in 
the coefficient functions fij characterizes the minimum 
photon energy and has been set to 5 = 0.9 [ pH , 
except for the first line in Eq. ( |2^ ) where 6—1.0 
corresponding to the full cross section. The f2i are 
complicated integrals which can be solved in terms of 
polylogarithms up to 5th order. In the code I use 
an expansion in 2: = m^/m^ and S — 1.0. Once 
experiments become more precise the correction to 
6 — 0.9 should be included. 

f{z) is the phase space factor for the semileptonic 
decay rate including NLO corrections [ I defined 
the MS mass ratio in z = [mc(mb)/mf,(r7?,b)]^ at the 
common scale, fi = rhh, which I also assumed for the 
factor TOj multiplying the decay widths. Since I do not 
re-expand the denominator this effects the phase space 
function at higher orders. Using the O(Q;3)-estimate0 

computed the 0(0^) coefficient for comparison only, and did 
not include it in the code. 
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from Ref. [ I obtain for the scmileptonic decay 
width, 



which I obtain the estimate. 



TsL ~ mlfo{z) 



1 + 2.7 ^ ' - 1.6 — 

TT V TT 



(23) 



where /o(z) is the leading order phase space factor. It 
is amusing that the coefficients in Eq. ( p^ ) are com- 
fortably (and perhaps somewhat fortuitously) small, 
with the 0(a3)-coeSicient even smaller than the one 
in Ref. [ where a low scale running mass had been 
advocated. Moreover, using the pre-factor in the 



numerator of R reduces the size of r-j in Eq. (|2^) and 
therefore the coefhcient k((5) — jTjiS) + rT/2 which 
multiplies the term as{Cj'^^Y . I obtain —2.1 < k{5) < 
1.4, while with the pole mass pre-factor one would 
have -8.7 < k{S) < -5.3. 
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A. Uncertainties from perturbative QCD 

Writing the perturbative expansion of some quantity 
in its general form for an arbitrary gauge group, it can 
easily be decomposed into separately gauge invariant 
parts. Table H shows for some (related) examples that 
after removing the group theoretical pre-factors, all 
coefficients, yi, are strictly of order unity, and that 
their mean, y, is very close to zero. In particular, 
there is no sign of factorial growth of coefficients. 
These observations offer a valuable tool to estimate 
the uncertainties associated with the truncation of the 
loop expansion, so I would like to make them more 
precise. 

Assume (for simplicity) that the yi are random draws 
from some normal distribution with unknown mean, 
/i, and variance, a^. One can show that the marginal 
distribution of fi follows a Student-t distribution with 
n — 1 degrees of freedom, centered about y, and 

with standard deviation. 



n{n — 3) 



(24) 



As can be seen from the Table, fi is consistent with zero 
in all cases, justifying the nullification of the unknown 
coeflticients from higher loops. I next assert that the 
distribution of cr, conditional on fi = 0, follows a scaled 
inverse-x^ distribution with n degrees of freedom, from 



a — ffQ ± Act 



i± 



2(n - 4) 



(25) 



Inspection of the Table shows indeed that cr, as the 
typical size of a coefficient, is estimated to be ^ 0{1). 

I now focus on the par tial hadronic Z decay width. 
As discussed in Section |2.3| , the 0(0;^) term, ^2, is 
much smaller than the tt^ term arising from analytical 
continuation. This is specifically true for the relevant 
case of n/ = 5 active flavors, where large cancellations 
occur between gluonic and fermionic loops. Notice, 
that the I?- function, in contrast to -Rhad, has opposite 
signs in the leading terms proportional to C\Cf and 
CaCfTfup. Indeed, the Adler D-function and the 
/3-function have similar structures regarding the signs 
and sizes of the various terms (see Table ||), and we 
do expect large cancellations in the /3-function. The 
reason is that it has to vanish identically in the case of 

= 4 supersymmetry. Ignoring scalar contributions 
this case can be mimicked by setting TpUf = 2Ca 
(there are 2 Dirac fermions in the = 4 gauge 
multiplet) or Uf — 12 for QCD, which is of the right 
order. In fact, all known QCD /3-function coefficients 
become very small for some value of between 6 and 
16. We therefore have a reason to expect that similar 
cancellations will reoccur in the di at higher orders. As 
a Icr error estimate for da, I suggest to use the largest 
known coefficient (3 x 0.71) times the largest group 
theoretical pre-factor in the next order {C'^Cf) which 
results in 

= ± 77. (26) 

With Eq. (|) and ds(Mz) = 0.120 one can absorb all 
higher order effects into the C'(ag)-coefficient of i?had, 
— —81 ± 63. This shifts the extracted as from 
the Z line-shape by -1-0.0005 and introduces the small 
uncertainty of ±0.0004. 

The argument given above does certainly not apply 
to the quenched case, n/ = 0, and indeed d2{nf — 0) is 
about —73% of the tt^ term, i.e., large and positive. In 
the case of = 3, which is of interest for the precision 
determination of from r decays, d2 is about —38% 
of the TT^ term. If one assumes that the same is true 
of 1^3, one would obtain d^(nf = 3) = 60. Estimates 
based on the principles of minimal sensitivity, PMS, 
or fastest apparent convergence, FAC, yield d'i{nf = 
3) = 27.5 [|2^ so there might be some indications for a 
positive d^{nf = 3). In any case, all these estimates lie 
within the uncertainty in Eq. ( |2^ ) and we will have to 
await the proper calculation of the C'(as)-coefficient to 
test these hypotheses. Note, that the current r decay 
analysis by the ALEPH Collaboration uses da = 50 ± 
50 [ pOl which is more optimistic. 
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The analogous error estimate for the five-loop (3- 
function coefficient yields, 



/34 = ± 579. 



(27) 



To get an estimate for the uncertainty in the RG 
running of ds , I translate Eq. (|27|) into 



P3 = P3 ± Pi, 



(28) 



where /ig is taken to be the lowest scale in- 
volved. This overestimates the uncertainty from /34, 
thereby compensating for other neglected terms of 
C'(a^+'' In" For the RG evolution from = 

to /ti = Mz this yields an uncertainty of Aas{Mz) = 
±0.0005. Conversely, for fixed a,{Mz) = 0.120, I 
obtain d^TOb) = 0.2313 ± 0.0006, d^(mr) = 0.3355 ± 

0. 0045. and as{mc) = 0.403±0.011, where I have used 
TOf, — 4.24 GeV and rhc — 1.31 GeV. For comparison, 
the ALEPH Collaboration quotes an evolution error 
of Aas{Mz) = ±0.0010 which is twice as large. I 
emphasize that it is important to adhere to consistent 
standards when errors are estimated. This is especially 
true in the context of a global analysis where the 
precisions of the observables enter as their relative 
weights. 
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Table 2 

Coefficients (MS ) appearing in the /9-function of a 
simple group [ ^; in non-Abelian corrections to the 
QED /3-function (denoted D) [ |, H; in the Adler 
I?- function [ ^ (rescaled by an overall factor 1/3); 
and in i?had (analytical continuation of D). The first 
four segments correspond, respectively, to the first four 
loop orders of non-singlet type. The fifth segment is 
the singlet (double triangle) contribution in 0{a'^). In 
D, D, and i?had, an overall factor aMz and the sums 
involving charges or Z couplings have been dropped. 
The completely symmetrical tensors of rank four, dA 
and dp, as well as Tp when appearing in parenthesis, 
apply to the /3-function only. Each Tp is understood to 
be multiplied by the number of fiavors, n/, except for 
the singlet term involving the symmetrical structure 
constants, d. 



group factor 


13 


D 


D 


-Rhad 


Ca 


0.92 








[-I-F) 


u. oo 




u.oo 


u.oo 


c'^ 

^A 


0.71 








CaTf 


-0.42 








Cf{Tp) 


-0.25 


-0.25 


0.25 


0.25 


^1 


0.83 








C\Tf 


-0.82 








CaCf{Tf) 
CaTI 


-0.36 


-0.23 


0.18 


0.18 


0.09 








CUTf) 


0.03 


0.03 


-0.03 


-0.03 


CfTf{Tf) 


0.08 


0.08 


-0.06 


-0.06 


C\ 


1.19 








C\Tf 


-1.67 








Cj^Cf{Tf) 

^A^F 

CaCp{Tf) 


-0.23 


-0.28 


0.32 


-0.38 


0.50 
-0.42 


-0.07 


0.51 


0.51 


CaCfTf(Tf) 
CaTI 


0.51 
0.01 


0.42 


-0.71 


-0.21 


CUTf) 

CITf{Tf) 

CfTUTf) 


0.18 


0.18 


-0.18 


-0.18 


-0.17 


-0.17 


0.02 


0.02 


0.02 


0.02 


0.09 


-0.01 


d^/NA 
dAdp/NA 


1.07 
-2.38 








dl/NA [Tld^A] 


0.50 


0.50 


-0.50 


-0.50 


y 


-0.02 


-0.01 


0.02 


-0.01 


A/i 


0.17 


0.09 


0.11 


0.09 


Co 


0.83 


0.28 


0.37 


0.31 


Ad 


0.13 


0.07 


0.09 


0.08 
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